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Abstract 

We construct jet observables for energetic top quarks that can be used to determine a short 
distance top quark mass from reconstruction in e + e~ collisions with accuracy better than Aqcd- 
Using a sequence of effective field theories we connect the production energy, mass, and top width 
scales, Q > m > T, for the top jet cross section, and derive a QCD factorization theorem for the top 
invariant mass spectrum. Our analysis accounts for: a s corrections from the production and mass 
scales, corrections due to constraints in defining invariant masses, non-perturbative corrections 
from the cross-talk between the jets, and a s corrections to the Breit-Wigner line-shape. This paper 
mainly focuses on deriving the factorization theorem for hemisphere invariant mass distributions 
and other event shapes in e + e~ collisions applicable at a future Linear Collider. We show that 
the invariant mass distribution is not a simple Breit-Wigner involving the top width. Even at 
leading order it is shifted and broadened by non-perturbative soft QCD effects. We predict that 
the invariant mass peak position increases linearly with Q/m due to these non-perturbative effects. 
They are encoded in terms of a universal soft function that also describes soft effects for massless 
dijet events. In a future paper we compute a s corrections to the jet invariant mass spectrum, 
including a summation of large logarithms between the scales Q, m and T. 
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I. INTRODUCTION 



Precise determinations of the top quark mass m are among the most important Standard 
Model measurements being carried out at the Tevatron, and being planned at the Large 
Hadron Collider (LHC) and a future International Linear Collider (ILC). A precise top mass 
determination is important for precision electroweak constraints, as well as extensions to the 
standard model like minimal supersymmetry pi]. The present combined measurement from 
the Tevatron is m = 171.4 ±2.1 GeV [0, [3] and mainly relies on methods where a number 
of top-mass-dependent kinematical quantities and observables are used in a global fit to 
determine the most likely top quark mass. For these fitting methods 0, |5| the observable 
most sentitive to the top quark mass is the top invariant distribution. It is obtained from 
reconstructing the total invariant mass of the top decay products. At the Tevatron the 
invariant mass distribution is being used in connection with other top mass dependent 
observables due to the limited statistics. 

In principle the reconstruction of the top invariant mass distribution provides the most 
natural way to measure the top quark mass since the peaked structure at the resonance is 
most closely related to the notion of the mass of a propagating massive and unstable degree 
of freedom. This method can be applied at the LHC and ILC where larger statistics are 
available. Experimental studies have concluded that at the LHC top mass measurements 
with uncertainties at the level of 1 GeV || [7]] can be achieved, while at the ILC even smaller 
uncertainties can be expected || [|. However, since the top quark is a parton carrying 
non-vanishing color charge, its mass is a priori not directly observable. In fact the top 
mass should be considered as a renormalization scheme-dependent coupling in the QCD 
Lagrangian rather than a physical object, just like the strong coupling a s . As such, the top 
mass obtained from reconstruction also depends on the method and prescription that is used 
to define the top invariant mass since the latter is not a unique physical quantity. In fact the 
notion of a physical particle whose squared four-momentum is the mass does not apply to 
the top quark if one asks for a precision in the mass value that is comparable to the hadronic 
scale. This is also reflected in a number of conceptual and experimental issues for top quark 
mass determinations that are associated with gluon radiation, underlying events, and the jet 
energy scales - effects that can never be fully eliminated for measurements of the top quark 
mass from reconstruction. Moreover certain top quark mass renormalization schemes are 
more suitable for precision measurements than others since the choice of scheme can affect 
the higher order behavior of the perturbative corrections as well as the organization of power 
corrections. Suitable quark mass schemes are compatible with the power counting and also 
lead to an optimal behavior of the perturbative expansion. Such schemes can be identified 
and defined unambiguously if the precise relation of the observable to a given Lagrangian 
top quark mass scheme can be established. 

For all jet based methods of top quark mass determination, and for reconstruction in 
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particular, these issues have been intrinsically difficult to address in the past. Previous work 
has not provided a coherent analytic framework in which perturbative and non-perturbative 
effects could be described in a systematic manner. Considering the expected precision for 
top quark mass measurements in the upcoming experiments such a framework is imperative. 

A top mass determination method where a systematic analytic framework exists and 
where the relation between the Lagrangian top mass parameter m and the measured top 
mass can be established to high precision is the threshold scan of the line-shape of the 
total hadronic cross section in the top-antitop threshold region, Q pa 2m, at a future Linear 
Collider [[HJ [11], [12], [13], |14|, where Q is the cm. energy. In this case the system of interest is a 
top-antitop quark pair in a color singlet state and the observable is related to a comparatively 
simple counting measurement. The line-shape of the cross section rises near a center of mass 
energy that is related to a toponium-like top-antitop bound state with a mass that can be 



computed perturbatively to very high precision [15], [U| [L7| [L8|, |19| using non-relativistic 
QCD (NRQCD) [5j| |l| an effective field theory (EFT) for nonrelativistic heavy quark 
pairs. The short lifetime of the top quark, r = 1/T pa (1.5 GeV) -1 , provides an infrared 
cutoff for all kinematic scales governing the top-antitop dynamics and leads to a strong 
power suppression of non-perturbative QCD effects. Experimental studies concluded that 
theoretical as well as experimental systematic uncertainties for this method are at a level 
of only 100 MeV p2], 531. The most suitable top quark mass schemes are the so-called 



threshold masses |Tj| , which can be related accurately to other short- distance mass schemes 
such as the running MS mass. Unfortunately, the threshold scan method cannot be used 
at the LHC because the top-antitop invariant mass can only be determined with a relative 
uncertainty of around 5% [EJj], which is not sufficient to resolve the top-antitop threshold 
region. 

In this work we use EFT's to provide, for the first time, an analytic framework that can 
be applied to systematically describe the perturbative and nonperturbative aspects of top 
quark invariant mass distributions obtained from reconstruction. As a first step towards 
developing a detailed framework for the LHC, we focus in this work on jets in a e + e~ 
Linear Collider environment at cm. energies far above threshold Q ~ 0.5 — 1 TeV. For 
e + e~ collisions strong interaction effects arising from the initial state can be neglected and 
there is no need to identify or remove any 'beam remnant' or underlying events. Also, in 
the e + e~ framework it is easier to formulate shape variables like thrust that control the 
jet-likeness and the soft dynamics of an event. We consider the double differential top and 
antitop invariant mass distribution, where each of the invariant masses, M t 2 and , are 
defined from all particles in each of the two hemispheres that are determined by the event's 
thrust axis. In Fig. [I] we show a sketch of such an event. Other invariant mass definitions, 
e.g. based on fc T algorithms and criteria to identify jets from top and antitop decay can be 
employed as well. Our approach also works for all-jet and lepton plus jet final states. Our 
focus is to study the double differential invariant mass distribution in the peak region close 
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FIG. 1: Six jet event initiated by a top quark pair, tt — * bWbW — * bqq'bqq'. The plane separating 
the two hemispheres is perpendicular to the thrust axis and intersects the thrust axis at the 
interaction point. The total invariant mass inside each hemisphere is measured. Our analysis 
applies equally well to the lepton+jets and the dilepton channels (not shown). 



to the top mass, so that M t 2 and Mf — m 2 ~ mT. It is convenient to introduce 

the shifted variables 

s t i= — = ~T<m, (1) 

m m 

because it is only the invariant mass distribution close to the peak that we wish to predict. 
Here the top width T is setting a lower bound on the width of the invariant mass distribution 
and the shifted variable s t) i can also be larger than T as long as s t j "C Tn. However, for 
simplicity we will often write s^j ~ T as we did in Eq. fll]). 

There are three relevant disparate scales governing the dynamics of the system, 

Q > m > T > A QCD . (2) 

This kinematic situation is characterized by energy deposits contained predominantly in 
two back-to-back regions of the detector with opening angles of order m/Q associated to 
the energetic jets coming from the top quark decay and collinear radiation. Frequently in 
this work we refer to the jets coming from the top and antitop quark collectively as top 
and antitop jet, respectively, but we stress that we do not require the jets from the top 
and antitop decay products to be unresolved as pictured in Fig. |1] (for example one can still 
identify a W and do 6-tagging). The region between the top jets is predominantly populated 
by soft particles with energies of order of the hadronic scale. 

The EFT setup used to describe the dynamics in this kinematic situation is illustrated in 
Fig. and represents a sequence of different EFT's. The use of different EFT's is mandatory 
to separate the various relevant physical fluctuations. The high energy dynamics for the 
top quarks at the scale Q 3> m can be described by quark and gluon degrees of freedom 
that are collinear to the top and antitop jet axes, and by soft degrees of freedom that 
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FIG. 2: Sequence of effective field theories used to compute the top/antitop invariant mass distri- 
bution in the peak region. 



can freely propagate between the jets. The appropriate EFT for this situation is the Soft 
Collinear Effective Theory (SCET) |25], F2(| |27|, with a nonzero top quark mass term 



which represents an expansion in A ~ m/Q ~ 0.2 — 0.3. The leading order soft-collinear 
decoupling |27j properties of SCET allows a factorization of the process into three sectors: 
top jet dynamics, antitop jet dynamics, and dynamics of the soft cross talk between the top 
and antitop jets, which corresponds quite intuitively to the situation pictured in Fig. |]. In 
SCET the typical fluctuation of the jet invariant masses around the top mass are still of 
order m, s t j ~ ni. Thus to describe invariant masses in the peak region s t j ~ T the top 
and antitop jets are finally computed in Heavy-Quark Effective Theory (HQET) [|30| which 
represents an expansion s/m and T/m ~ 0.01. We have in fact two copies of HQET, one 
for the top and one for the antitop, plus soft interactions between them. In these EFT's the 
top decay can be treated as inclusive and is therefore described by the total top width term 
T that acts as an imaginary residual mass term [T(], |3TJ] . Since HQET is usually understood 
as being formulated close to the rest frame of the heavy quark without the soft cross-talk 
interactions, we refer to these two EFT's as boosted HQET's (bHQET's). 1 

At leading order in the expansion in m/Q and T/m we show that the double differential 
invariant hemisphere mass distribution can be factorized in the form 

( rfM^M- 2 ) h = a o H Q(Q^Q^rn)H m (m,Q,fi m ,^ (3) 

x / d£ + d£-B + (s t - T, fj) B_ (s- t - 9L } r, /A S hemi (l+, £~, //) , 
J \ m J \ m J 



We adopt the acronym bHQET in cases where we wish to remind the reader that the residual momentum 
components of the heavy quark in the e + e~ cm. frame are not homogeneous, and that additional gluon 
interactions occur which are not simply the soft gluons of standard HQET. 
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where §t and % are denned in terms of M t 2 f in Eq. (JJ). The term cq is a normalization factor, 
and the factors Hq and H m are matching corrections that are derived from matching and 
running in SCET and the bHQET's, respectively. Hq and H m are independent of s t and % 
and do not affect the form of the invariant mass distributions. The jet functions B± describe 
the QCD dynamics of collinear radiation in the top/antitop direction, and the decay of the 
top and antitop quarks near mass shell within the top/antitop jets. They can be computed 
perturbatively at the scale fi > T since the top width T provides an infrared cutoff from 
hadronization. At tree level they are Breit-Wigner functions 

B ± (s,T) = + ... , (4) 

Tim + L 1 

where the ellipses indicate QCD corrections that distort the Breit-Wigners. For the compu- 
tation of the B± it is mandatory to employ properly defined short- distance top mass schemes, 
to achieve a well-behaved perturbative expansion. Finally, the soft function Sh em i(£ + > &~) 
describes the physics of the soft nonperturbative gluons through which the top and antitop 
jets can communicate. The low energy fluctuations of these soft gluons are not cut off by 
the large top quark width. This can be intuitively understood due to the lifetime dilation 
of the boosted top quarks. As explained in Sec. |I[ using soft-collinear factorization, we 
can show that the soft function is universal, namely that the same function governs the low 



energy dynamics for massless jets in the dijet limit [[!2|, [33], [34], |35], |3"6fl . So, information on 
the form of Shemi(^ + , can be gained in a model-independent way from experimental data 
on massless dijet events. The form of the factorization theorem in Eq. (13) is based on the 



same principles as the factorization formula for massless dijet event shapes |32], |33|, |34], |36 
but it differs due to the need to treat massive quark jets and effects related to the large top 
quark width. We also use our results to derive a factorization theorem for thrust and the 
heavy-jet mass event shape for ti production in the peak region. These distributions can 
also be used to measure the top-quark mass. 

The convolution in Eq. (|||) shows that the observed hemisphere mass distributions are 
inevitably distorted by the nonperturbative soft momentum distribution, and that top and 
antitop jets can only interact indirectly through exchange of different light-cone momentum 
components that are governed by the soft function. We can also show that for invariant 
masses M t) j that are defined through the identification of the jets from top and antitop 
decay, which are determined from a kx jet algorithm, the same factorization formula as in 
Eq. ([D can be derived up to a different soft function. 

We believe that the factorization approach proposed in this work, and the factorization 
formula in Eq. (|3|), represent advancements concerning the following points: 

• We give a well defined relation between a jet observable sensitive to the top mass and 
the Lagrangian mass. This allows the definition of a short- distance top mass which 
we call the "jet-mass" . Theoretically the jet-mass can be determined with a precision 
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better than Aqcd, once the soft function governing nonperturbative effects is known 
by other means. We expect that the jet-mass will be useful for a broad range of 
observables involving jets and parton showering from massive quarks. 

• The soft function appearing in the massive-jet factorization formula is universal, and 
appears in massless dijet event shapes. This universality can reduce the dependence 
of uncertainties in the top mass from reconstruction on parton shower monte carlos 
and hadronization models. 

• The factorization approach opens up the possibility to systematically construct top 
mass observables where nonperturbative effects are suppressed. 

While the focus of this paper is on tt production at a e + e~ Linear Collider, the main ideas 
and tools developed are general, and will also play an important role for the environment 
of the LHC where a substantial number of top events with large pt will be available. Other 
applications of our approach to factorization of jets from massive particles may include 
processes such as single top production |37], [38J, W pair production, or processes involving 
new colored unstable particles |39|, |3(J. We briefly comment on these applications in the 
summary. 

The outline of the paper is as follows. In Sec. [Tl] we describe the relevant EFT formalism 
for our computation. In Sec. |T| we derive the factorization theorem in SCET, introduce the 
hemisphere jet invariant masses, and perform the factorization of mass effects in boosted 
HQET. The result of the analysis in this section is the complete factorization theorem for 
the double invariant mass distribution, and the extension to thrust and the heavy-jet mass 
event shapes. In this section we also define the short- distance jet-mass scheme. In Sec. [TV] 
we study the factorization theorem numerically at leading order, and discuss implications for 
top-mass measurements. We also display numerical results for the shape of the peak region. 
In Sec. [V| we discuss the relation between the factorization theorem for the hemisphere 
invariant masses used in our work, to a factorization theorem for the reconstruction method 
based on kj< jet algorithms employed in Refs. 0, M. Finally we summarize and conclude in 



Sec. VI. 



This paper concentrates on the derivation of the factorization theorem, on field theoretic 
issues, and on the basic phenomenological implications of our result. Readers only interested 
in the final result may skip over the analysis in sections [TI| and |TTTJ, and go directly to 
section [TV]. In a future paper we present the computation of a s corrections to the jet invariant 
mass cross-section, and the summation of large logarithms between the scales Q, m, T. 
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II. THE EFFECTIVE FIELD THEORIES 



In this section we discuss the EFT's required to compute the double differential invariant 
mass distribution d 2 a/dM 2 dM 2 in the peak region. The relevant energy scales are: 

Q > m > s t ~ % ~ T , (5) 

where the hatted s- variables were defined in terms of M 2 and M| in Eq. (|l|). Once radiative 
corrections are included, large logarithms arise through ratios of the above energy scales, 
some of which are double logs, and thus can be quite large. For example Y/m ~ 1/120, so 
ln 2 (r/m) « 25. It is obviously important to understand the appearance of all large logs as 
accurately as possible, and to sum them systematically. This summation is accomplished 
by matching onto a sequence of EFTs and using renormalization group equations (RGE's). 

Starting from QCD we first switch to the Soft Collinear Effective Theory (SCET) |2!| 
27 , |28| , 4lj for massive quarks and then to Heavy Quark Effective Theory (HQET) ||42| , |43| , 



44 , |i"5| , |46|| combined with the unstable particle EFT method |3l| , [47| , |48| , j4"9"f . This scheme 
includes systematically effects related to the large top quark width, as well as interactions 
related to the soft cross-talk: 

QCD — ► SCET — ► boosted-HQET with unstable heavy quarks. (6) 

An intuitive picture which displays why this sequence of EFTs is relevant is shown in Fig. |2|. 
We are interested in events where the top quarks are produced close to their mass shell as 
characterized by the condition in Eq. (p]). At the production scale Q, the invariant mass 
of the top and antitop quarks can still fluctuate with s t) j ~ Q due to it's interactions with 
hard gluons of characteristic momentum ~ Q. In the first step, when switching to SCET, 
these hard modes are integrated out and we expand mm/Q -C 1. SCET makes it simple to 
separate the physics associated with i) the top-jet, ii) the antitop jet, and iii) the soft-cross 
talk between the jets. After the implementation of this factorization theorem, each top jet 
and the soft-cross talk can be studied independently in the field theory. The factorization 
theorem tells us how to tie them together. Now in SCET the invariant mass of the top quark 
fluctuates with s t) t ~ m , so we still have to remove these large momentum fluctuations to 
describe the desired kinematic region where Y ~ s t f <C m. Such invariant mass fluctuations 
are analogous to those encountered in HQET for a bottom quark inside a 5-meson 

(mv + k) 2 — m 2 = 2mv ■ k + k 2 ~ 2mAqcD , (7) 

with the difference that for the unstable top quark v-k — > v-k+iY/2. Since top-quarks decay 
before they have a chance to hadronize, the top-width Y adopts the role Aqcd plays for the 
B meson. Keeping in mind that the tops are highly boosted and unstable, we actually match 
onto two boosted versions of HQET, one for the top and one for the antitop. A discussion of 
the necessary SCET and HQET theoretical ingredients is given in the following subsections. 
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bCEl [A ~ m/Q <C 1J 


briQEl [1 /m <C 1J 


n-collinear (£ n , A%) p^~(2(A 2 ,l,A) 
n-collinear (&>, A£) ^~Q(1,A 2 ,A) 
mass-modes (q m , A&) p^ ~ Q(X, A, A) 
Crosstalk: soft (q s , AS) K~Q(A 2 , A 2 , A 2 ) 


n-ucollinear (/i„ + , A£) A:' 1 ~T(A, A -1 , 1) 
n-ucollinear (h v _ , At) ~r(A~ 1 , A, 1) 

same soft (q s , As) ps ~(A,A, A) 



TABLE I: Summary of the fields required in SCET and bHQET. The first field in each bracket is 
a quark, and the second is a gluon. The scaling of momentum components is given for (p + ,p _ ,p _L ). 
After factorization, the soft fields on the last line generate a cross-talk theory that communicates 
with collinear fields in both SCET and bHQET through two kinematic variables. A is the scale 
for the soft modes. 



A. SCET with Masses 

SCET is an effective theory describing the interactions of soft and collinear particles, 
which are characterized by the scaling of their momenta. In this framework it is convenient 
to introduce the four-vectors 

n^ = (l,n), r^ = (l,-n), (8) 

where n can be thought of as the direction of the top jet and — n as the direction of the 
antitop jet (n 2 = 1, n 2 = 0, n 2 = 0) . Any momentum can then be decomposed as 

P* 1 = n-p — + n-p — + p^, (9) 

and we denote momentum components in this light cone basis as {p + ,p~ ,Pj_) = {n-p, n-p, p±). 
The square of the momentum vector p M then reads p 2 = p + p~ +p\- It is also convenient to 
denote the momentum of collinear particles in the n and — n directions by the subscripts n 
and n respectively, which corresponds to the large energy modes in the corresponding jets. 
Thus we have collinear labels 

n for the top-jet, n for the antitop-jet . (10) 

The momentum of soft particles that communicate between the jets will be denoted by a 
subscript s. We also have mass-modes that are required in order to describe certain top- 
quark vacuum polarization loops. The momenta of the collinear, mass, and soft modes 2 have 
the typical scalings shown in table |in the SCET column, where here A is the small expansion 
parameter. A particle with components scaling as (A 2 , 1, A) has a small A-momentum relative 



2 In some factorization theorems it is necessary to distinguish between soft and ultrasoft particles, and 
between two versions of SCET: called SCETj and SCET n - In this paper we only deal with SCETj with 
ultrasoft gluons. For simplicity we will therefore simply use the term soft modes. For modes with momenta 
~ (to, to, to) that are specific to the massive SCET theory, we use the term "mass-modes" . 



10 



to its energy, and is said to be collinear to the n M direction etc. Both A and the hard scale Q 
have a size that depends on the particular process under study. For example, in B — > X s 7 
the hard scale is the 6-quark mass mj, and the expansion parameter is a/AqodT^- For 
pair production of top jets, the hard scale Q is the center of mass energy, and the SCET 
expansion parameter is 

A~|. (11) 
It follows that the typical virtuality of the collinear, mass, and soft modes in SCET satisfy 



Pl~P 2 n~ ™ 2 , P 2 m ~ m2 i and P 2 ~ ( 12 ) 

Since m 4 /Q 2 3> Aq CD , the soft modes in this theory still contain perturbative components 
as well as the underlying non-perturbative dynamics at smaller scales. Using m = 171 GeV 
this is true for Q <j 40 TeV i.e. for any conceivable cm. energy of a future Linear Collider. 
The soft particles correspond to modes with wavelengths that allow cross talk between the 
two jets. In addition, at two-loop order the soft gluons in SCET interact with virtual 
top-quarks which are described by the mass-modes indicated in table |. These mass-modes 
do not interact directly with the collinear fields and only appear as virtual effects for our 
observable, because we only consider cases where s t j <C Qm. (As discussed below Eq. (p3|) .) 
In addition we have virtual collinear top-quarks that can interact with the collinear particles 
through the collinear Lagrangian. The n-collinear, n-collinear, mass modes, and soft modes 
are described by separate quark and gluon fields which are also listed in table |. Hard modes 
involving momenta p" ~ Q have already been integrated out when QCD is matched onto 
SCET. 

At leading order the SCET Lagrangian for collinear particles in different directions can 
be written as a soft Lagrangian plus a sum of collinear terms HTJ, = C s + Cn}- 



i 



The sum satisfies the constraint • rij ^> A 2 for i j£ j, with the choice of A determining 
what is meant by distinct collinear directions. The collinear particles in different sectors 
only interact via soft gluon exchange or interactions in external operators. When the _L- 
momentum of the collinear particles is of the same size as the quark mass the result for the 
leading order collinear Lagrangian [^6, ^7]] must include the quark mass terms derived in 



Ref. (see also Ref. ]50]). The collinear quark Lagrangian for the direction n is therefore 
given by 

45 = [in ■D s + gn-A n + (i^-m)W n ^W^+m)] t( n , (13) 



with D^r ~ m 3> D^. There is also an n-collinear Lagrangian for gluons [28]. Here the soft 



and collinear covariant derivatives are 

iD* = id" + gA^, iD£ = V" + gA% , (14) 
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soft particles 




mass is restricted and the top-decay products become explicit by matching onto HQET. 



where V 1 is a label operator picking out the large collinear momentum of order Q and QX of 



a collinear field |27| , while the partial derivative acts on the residual momentum components 
<9 M ~ A 2 . The term W n is the momentum space Wilson line built out of collinear gluon fields 



W n {x) = ^ exp ( - • A n {x) ) . (15) 

perms 

We also note that Eq. flT5| ) is the bare Lagrangian. In particular, any mass definition can 
be chosen for m through an appropriate renormalization condition without breaking the 
power-counting. At 0(a s ) these mass-schemes are the same as those in QCD [pT]| , because 
the self-energy graphs are directly related. 

An example of an external operator that connects different collinear sectors is the jet 
production current, which couples to the 7* or Z*. In QCD the production matrix element 
is (Xl^jO) where (X\ is the final state. The required vector and axial currents are given 
by 

J*(x) = $(x)W(x) , J?(x) = ^(x) 7 M 75^) , (16) 

and for convenience we will adopt the short-hand notation = ip(x)r^ip(x). The matching 



relation of these QCD currents to SCET currents is given by the convolution formula [26] 



Jf(0)= [duduC{u,u,{M)4 0)lx (u,u>,(t), (17) 

where C contains short- distance dynamics at the scale Q, while jf^* 1 describes fluctuations at 
all longer distance scales. In the presence of multiple collinear fields, as well as modes scaling 
like our mass-modes and soft-modes, the construction of currents in SCET has been discussed 
in great detail in Ref. JO]. Interactions between the mass- modes and the collinear- modes 
produce offshell particles, which when integrated out leave residual interactions through 
Wilson lines in the SCET current. The SCET production current at leading order in A is 
given by 

4 0) >, O, fi) = xnA^sirfs n xnAo) , (is) 
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where X nA Q ) = % ~ n-V)(W^ n )(0) and X nd Q ) = S(p ~ n-V)(W^ n )(0). The mass- 
mode Wilson lines and S n will be described below. Here the (0) indicates that the fields 
are at coordinate x M = 0, and we recall that this x M dependence carries information about 
the residual momenta at the scale QX 2 = m 2 /Q. The dependence on larger momenta is 
encoded in labels on the collinear fields |27[], and, for example, S(u — n ■ P) forces the total 



minus-label-momentum of (W^ n ) to be uj. We also use the notation Xn = (W^ n ) and 
One can decouple the soft and collinear modes in Cq°2 by performing a field redefinition 



on collinear fields 

£n * Y n ^ n j * Y n J (^) 

where F n is a soft Wilson line 

Y n (x) = P exp ig J dsn-A s (ns + x)^j . (20) 

This gives 

Y%(x) = P exp dsn-A s (ns+x)^ , (21) 



which satisfies Y^Y n = 1. For two-jet production the factorization is most transparent |34[ 



with the reference point sq = oo shown in Eq. ( ^0|) . The gluon fields are either antipath- 
ordered (for P) or path-ordered (for P). We use the same Wilson line for both the quark 
and antiquark parts of £ n . Another possibility is to make different field redefinitions on the 
particle and antiparticle parts of the fields [|52]| . In fact, all results are independent of the 
choice of reference point in the field redefinition; the path is determined entirely by changes 



the field redefinition induces on the operators and the interpolating fields for the states ||53 

The mass-mode Wilson line S n (x) is defined in an identical manner to Eq. (|2ll|), but 
with n ■ A s — ► n ■ A m . In order to avoid double counting with the effects contained in the 
soft Wilson lines the mass-mode Af^ fields are defined with zero-bin subtractions for the 



soft-region p>l, and we have mass-mode top-quarks ip m with a mass m. Any graphs with 



mass-mode gluons that do not involve a top-bubble from ip m fields are exactly canceled 
by these zero-bin subtractions. Thus the mass-modes only contribute in these vacuum 
polarization graphs. The soft-gluons can also couple to the ij) m fields, however they do so 
with a multipole expansion, and therefore do not inject momentum into the closed ip m loop. 

After the change of variable in Eq. (|19D the leading order SCET collinear quark La- 
grangian and current become 



4? =L[in-d + gn-A n + {i^-m)W n iw^(i^+m) 



nq 
r(°)M 



2 Zn 



Ji ^ — Xn,uXlSl^^S n Y n Xn,u){ty > (22) 
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where we used the property Y^n ■ D s Y n = n ■ d. The only coupling of the soft gluon to the 
collinear quark was through in ■ D s which is no longer present (and a similar property occurs 
in the collinear gluon action). These soft couplings reappear as Wilson lines in the current 
as shown above. Hence we have achieved soft-collinear decoupling in the Lagrangian and 
the current. 

In the two-jet process we wish to factorize we must also consider the transformation 
property of the state |X) under Eq. fllpp. For manipulations in the factorization theorem 
for two-jet production we decompose the state into collinear and soft pieces, 

(X | = (X n X n X s \ . (23) 

Note that this decomposition is only valid for the states we are interested in for describing the 
dijet region, not for a general state in QCD. Since there is always at least one n-collinear and 
one n-collinear particle, we do not consider any mass-modes, X m , in these states either. The 
presence of a mass mode would induce an invariant mass p 2 x = (p n + Pm) 2 — VnPm — Q m ^ 
m 2 , which would make it impossible to satisfy the invariant mass condition required to study 
the peak region. Therefore the mass-modes will only appear as virtual contributions. The 
collinear states (X n \ and {X n \ are a color triplet and color antitriplet, just like a quark and 
antiquark state. Therefore, we must consider how these collinear states transform under 
the change in the action induced by Eq. (0). However, because these color triplet states 
can be derived from the out states at large time, t — » oo, they are not affected by the field 
redefinition with reference point at oo 5"3fl . With the current at x, we therefore have 



(X\T{x n ^SirS nX n,*}\0) - (X\T{x n ^StrS n Y n Xn,*}\0). (24) 

Here the T reminds us to keep the proper time-ordering of the A a (x) gluon fields in the Y's. 
There is no ordering issue between fields in Y n with those in Y^, since they are space-like 
separated and commute [M\. We also need the complex conjugate of Eq. ([24]) for the matrix 
element, which is 

(0\T{x n , a SlrS n Xn,.}\X) -> (0|T{x^ StT S n Y nXn ,u}\X) , (25) 
where T is anti-time-ordering. Note that 

/ /'OO . 

T (Y fl ) T = (Yi)* = % ] = P exp [ig J ds n-A s (ns+x)) , 

T (Y^) T = (Y n y =%=P exp(-igj dsn-A s {ns+x)^) , (26) 

where A s = AfT A , with T A = —(T A ) T the generator for the 3 representation, and the super- 
script T is the transpose with respect to the color indices of the fundamental representation. 
If we switch to these barred Wilson lines then the time-ordering and anti-time-ordering be- 
comes redundant. Eq. (|26|) applies equally well for the S Wilson lines. Considering the 
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squared matrix element for the cross-section we find 

(o\T {xn,^ strs n Y nXn ^}\x)(x\T{x n ^ si r s n f^}|o) 

= (0\TnA Y *) ba (Sn) b ' b (TS n Y nXn ^) b '\X)(X\(x n ^ ^r) c ' (ft)* (f*)*^) 

= Mfari (o\TnA y n) ba (rY nXn ^) b '\x)(x\(x n ^ry' (Ft) dc 4Jo> , (27) 



where a, b, c, d, d are color indices. The decoupling of soft gluons in Eq. (p7|) is identical to 
that in massless two-jet production, and ignoring the mass-mode Wilson lines the discussion 



above agrees with the SCET derivation in Ref. [53], as well as the original derivation in 



Refs. |)5, 56]. To obtain the last line in Eq. (|2"7|) we note that the Dirac structures V and 
r are color singlets, and that the mass-mode Wilson lines can be separated into vacuum 
matrix elements since there are no mass-modes in the states. Furthermore 

(0\(S n ) b ' b (S n ) b ' a \0) = 5 ^(0\(S n ) bW (S n ) b ' a '\0) , (28) 

with an analogous result for (0|(S'^) ac '(S'^) cc '|0), so this contracts the color indices on either 
side of the product of soft Wilson-line factors. Thus defining 

M{m,ri = ±\{0\-3ZS*\0)\ 3 , (29) 

we are left with the matrix element shown on the last line of Eq. fl27|) . Here A4(m,fi) = 
1 + 0(0?). 

The soft-collinear decoupling property is crucial to organizing the physics of the massive 
two-jet problem. As we will show in Sec. [TTI|, there is a factorization theorem in SCET that 
decouples the soft and collinear modes at leading order which allows us to study the physics 
of each jet independently. The cross-talk is confined to a simple top mass independent 
vacuum matrix element involving the Y Wilson lines, 

(0\(Y n ) ba (Y n ) bc \X s )(X s \(Y^) cd (Yl) ad \0) , (30) 

which agrees with the corresponding soft matrix element for massless quark production |36|, 



56], |5jJ and which will eventually determine the soft function S , h em i(^ + , to be used in 
Eq. (|). As we also show in Sec. |T|, the precise definition of the soft function S depends on 
the prescription that is used how the momenta of soft particles enter the top and antitop 
invariant masses §t and %, respectively. In the next subsection we describe how the matrix 
element in Eq. ([30]) is modified when we integrate out the top-quark mass. 

Finally, in SCET because the top-quark mass m, and mass of the H^-boson, m^y, are 
still low energy scales, the decay of an n-collinear top-quark is simply described by the full 
electroweak interaction, 

C ew = ^= bW->fP L t + ^= tW+^PLb , (31) 
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where Gp = \[2g\/{^m 2 v ) is the Fermi constant. This treatment is consistent since we can 
treat the top decay as fully inclusive up to 0{m 2 /Q 2 ). 

The collinear A n and A n gluons in SCET can induce fluctuations s t , % ~ m. Once we 
restrict ourselves to events with T ~ s t , sj <C m, i.e. we force the top quark and antiquark 
to remain close to their mass shell, the situation looks very much like two distinct copies 
of HQET in boosted frames. There was nothing special in the dynamics that sets the scale 
m 2 /Q for the soft interactions, and so we call A the scale that controls the soft-cross talk. 
In the field theory A will be defined as the scale where we model or fit the primordial 
soft function. Generally we will take m ^> A ~ T, although any value A > Aqcd can be 
considered. So we must switch from SCET onto these HQET theories, and also consider 
what happens to the decay interaction in Eq. (|3lD . We describe the boosted HQET theories 
in detail in the next section, and we also discuss how the soft cross-talk interactions remains 
active when the fluctuations at the top mass scale m are integrated out. 

Since the above Lagrangians and currents are LO in A, it is natural to ask about the 
role of power corrections. As it turns out, higher order Lagrangians and currents give 
corrections to our analysis at 0(a s m/Q), 0(A/Q), 0{m 2 /Q 2 ), or 0{T/m). The absence 
of 0(m/Q) implies that the m/Q expansion does not significantly modify the top-mass 
determination. The leading action contains all m/Q corrections that do not involve an 
additional perturbative gluon, so the corrections are 0(ce s m/Q). We have also verified 
that at tree level the m/Q corrections to the SCET current [^DJ vanish when contracted 
with the leptonic tensor. Furthermore, many of the higher order m/Q corrections have the 
form of normalization corrections, and thus do not change the shape of the invariant mass 
distribution. Subleading soft interactions are 0(A/Q). The interplay of our hemisphere 
invariant mass variable with the top decay can induce 0(m 2 /Q 2 ) corrections, as we discuss 
later on. Finally there will be power corrections of 0(T/m) in bHQET. 



B. Boosted HQET with Unstable Particles and Soft Cross- Talk 



Boosted Heavy Quarks. HQET B2L 43, M, 4b| is an effective theory describing the 



interactions of a heavy quark with soft degrees of freedom, and also plays a crucial role for 
jets initiated by massive unstable particles in the peak regions close to the heavy particles 
mass shell. The momentum of a heavy quark interacting with soft degrees of freedom can 
be written as 

p» = m / + (32) 

where denotes momentum fluctuations due to interactions with the soft degrees of freedom 
and is much smaller than the heavy quark mass \k^\ Cm. Also typically v ^ ~ 1 so that we 
are parametrically close to the top quark quark rest-frame, = (1,0). 
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In the top-quark rest frame we have r ^ <C m, where refers to momentum fluc- 
tuations of the top due to interactions with gluons collinear to its direction, which preserve 
the invariant mass conditions T ~ §t, % -C m. For our top-quark analysis, the center of 
mass frame is the most convenient to setup the degrees of freedom. In this frame the gluons 
collinear to the top-quark which preserve the invariant mass condition will be called ultra- 
collinear (ucollinear) in the n direction. A different set of n-ucollinear gluons interact with 
the antitop quark which moves in the n direction. The leading order Lagrangian of the EFT 
describing the evolution and decay of the top or antitop close to it's mass shell is given by 

— % — % 

C + = h v+ (iv + ■ D + — 5m + -T) h v+ , £_ = h v _ {iv _ ■ L>_ — 5m + -T) h v _ , (33) 

where the + and — subscripts refer to the top and antitop sectors respectively, and iD^ = 
i<9 M + gA±. These HQETs represent an expansion in T/m. The HQET field h v+ annihilates 
top quarks, while h v _ creates antitop quarks. In the cm. frame the components of fc M are 
no longer homogeneous in size, and t>± Y 1. Instead for the (+, — , _L) components we have 



it 

v_ 



m Q \ , ,, ( m Q , 
Q m J \Q m 1 

Q m \ , ,, ( Q m . 

-,^,0, ), F ~ r( — , — , i ]. 



m Q J \m Q 



Note that the T in Eq. ( p4|) can be replaced by a larger scale, of order s, as long as this 
scale is much less than m. Eq. (0) is easily obtained by boosting from the rest frame of 
the top and antitop respectively with a boost factor of Q/m. In this naming scheme we will 
continue to call the gluons that govern the cross-talk between top and antitop jets soft. We 
emphasize that they are not included in £±, since they have nothing to do with the gluons 
in standard HQET. Soft gluon interactions will be added below. To avoid double counting 



between the soft gluons, the ultracollinear gluons are defined with zero-bin subtractions |54 
so that for example n-k + ^ and ^ 0. Finally, since HQET is applied for fi < m there 
are no analogs of the SCET mass-modes in this theory. All effects associated with virtual 
top-quark loops are integrated out at the scale m. 

The leading order Lagrangians C± contain a residual mass term 5m which has to be chosen 
according to the desired top quark mass scheme. For a given top mass scheme m, the residual 
mass term is determined by its relation to the pole mass m po i e = m + 5m. Anticipating 
that we have to switch to a properly defined short- distance mass definition [J5S], ^ |50|, |5T[ ] 
when higher order QCD corrections are included, we note that only short- distance mass 
definitions are allowed which do not violate the power counting of the bHQET theories, 
5m ~ T. This excludes for example the use of the well known MS mass, since in this scheme 
5m ~ a s m ^> T. In practice, this means that using the MS mass leads to an inconsistent 



perturbative expansion as explained in Sec. [111 G| . This is the reason why the MS mass can 
not be measured directly from reconstruction. 
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The leading order Lagrangians C± also contain top- width terms iT /2. An effective field 
theory treatment of the evolution and decay of a massive unstable particle close to its mass 
shell was developed in |12|, [19], |31], (47], |lpj, fE| . The examples treated in these references were 



the resonant production of a single unstable scalar particle, and the leading and subleading 
width corrections to threshold ti production. In our case, we deal with the energetic pair 
production of massive unstable fermions, and we arrive at two copies of this unstable HQET 
corresponding to the top and antitop sectors. In these two HQET theories we treat the 
top and antitop decays as totally inclusive, since we do not require detailed differential 
information on the decay products. So the total top width Y appears as an imaginary mass 
term in £±, which is obtained by simply matching the imaginary part of the top and antitop 
self-energy graphs from SCET onto bHQET. As we show in Sec. p| , this inclusive treatment 
of the top decay is consistent with the hemisphere invariant mass definition we employ in 
this work up to power corrections of order (m/Q) 2 . We will come back to the role of higher 
order power corrections in the treatment of the finite top lifetime at the end of this section. 

Soft Interactions. Lets consider how the soft gluons interact with our heavy quarks in 
each bHQET. For a heavy quark in the boosted frame we consider interactions with soft 
gluons of momentum 

F<(A,A,A). (35) 

Our main interest is in the case A < T, but it is useful to keep a more general A > A QCD for 
the moment. We wish to demonstrate that these gluons are still entirely described by the 
cross-talk matrix element in Eq. (|30|), and that this is true without needing to expand in the 
ratio of A to T. Or in other words, that the simple eikonal propagators for the soft-gluon 
attachments to the energetic tops remains valid even below the mass of the quarks and even 
in the presence of the top-width. Our demonstration assumes the reader is quite familar 
with Ref. [28]. To prove this we go back to the original SCET Lagrangian in Eq. flT3] ) prior 



to the field redefinition, and match the soft interactions onto the HQET theory. This gives 
the same Lagrangian as in Eq. (^) but with replacements 

iD » _> iV % = id* + gA% + —gn- A s , 

iDt -> = idt + gAt + — gn- A s . (36) 

The new covariant derivatives T>± also appear in the pure gluon action responsible for the 
ultracollinear gluon kinetic term. The nature of the expansion for different momenta in id 11 
will depend on the size of the soft scale A relative to the smallest ultracollinear components 



mY/Q displayed in Eq. fl34j) . Note in this comparison that the width is suppressed by a 
factor of m/Q. Physically this factor is easy to understand, it is simply the time-dilation of 
the width of the energetic top-quark from the point of view of the soft gluons. The boost 
factor is encoded in v+ ■ n = t>_ • n = Q/m. 
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For our analysis we also need the effective current in the bHQET theories that corresponds 
to the SCET current in Eq. It is 

J bHQET = (h v+ W n )Vl;(wlh v _) , (37) 

where the Wilson lines are the same as W n and in SCET, except here we have gluons 
n- A + with path along n M for W n , and n-A_ with path along n M for W^. The simplest way to 
derive this result is to note that the two collinear sectors in the SCET current in Eq. fll8|) do 
not directly interact, and neither do the two sectors of the two bHQETs. In the rest frame of 
the top-quark for example the matching is simply (W%if) s ) (W^h v+ ), where ip s is a field for 
the top-quark near its rest-frame. Boosting this result gives the matching for the top-quark 
field in Eq. (|37"D, and the result for the antitop quark is analogous. The dynamics of the B + 
and B_ jet functions will be defined by the two interpolating field operators (W^h v+ ) and 
(W^h v _), and is governed by the Lagrangians C + and £_ respectively. 

Let us now come back to the derivation of Eq. (|30|) from bHQET. For convenience we 
start by taking both scales the same size, mT/Q ~ A. (Below we will show that the same 
result is obtained for the case where mT/Q <C A, which includes the situation A ~ T.) For 
mT ~ QA we can formulate the multipole expansion for the coupling of soft gluons to the 
heavy quarks by splitting the momenta into a large label components of size QT/m and T, 
and residual momentum components of size mT/Q. Thus 3 

z3£ = —n-V c + + —n-id , i&t = —n-V c + + —n-id . (38) 

The notation indicates that soft momenta only appear in the components id 11 . On the other 
hand ultracollinear momenta appear in all four components, and are picked out by the label 
operators V[? or the id. Next we make the same field redefinition on bHQET fields that we 
made on the SCET fields in Eq. (|T9|) 

h v+ > Y n h v+ , > Y n A^_Y^ , h v _ ► Yfji v _ , A^L > Y n A^_Y n , (39) 

where the fields in Y n and Y n are soft gluons. Since 

(v+-n)(in-d+gn-A a )Y n = Q, (v_-n)(in-d+gn-A s )Y n = 0, (40) 

this field redefinition gives back exactly Eq. ([£J) for the bHQET Lagrangian and also gives 
a leading ucollinear gluon action that has no couplings to soft gluons. In addition when 
making the field redefinition in the bHQET currents, Eq. fl37|) , we get exactly the same soft 
cross-talk matrix element for the two-jet production 

(o|(F,) 6a (^) bc |x s )(x s |(y;t) cd (Ft) ad |o). (4i) 



3 This formulation of the multipole expansion is the same as for the coupling of ultrasoft particles to 
collinear particles in SCET |28| where the two types of derivatives are formally separated by introducing 
label operators, and leaving residual momenta to be picked out by id 11 . 
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The only difference between the SCET matrix element in Eq. ( p0|) and the HQET matrix 
element in Eq. ([II]) is that in the former the soft-gluons couple to the massive ip m fields, while 
there are no such couplings in the latter. In matching renormalized soft matrix elements at 
a scale fi ~ m the only effect of these couplings to ip m fields is to induce an overall Wilson 
coefficient, so that S SCET = T (m, /i)S' bH Q ET . Thus the main dynamics of the soft gluons is 
not modified in a substantial way by passing from SCET to the boosted HQET Lagrangian, 
nor by the presence of the width term for unstable quarks. 

For completeness, let us now consider the case A ^> mT/Q and show that the same 
result is obtained. In this case a soft gluon of momentum coupling to an h v+ with 
residual momentum k + , has n ■ £ ^> n ■ k + , while for h v _ we have n ■ i ^> n ■ k-. Thus these 
soft gluons knock the heavy quarks away from their mass shell, and their interactions can not 
be formulated in a local manner in the same theory as the ucollinear gluons. This is similar 



to how soft and collinear gluons interact in the theory SCETn as discussed in Ref. [28 



To derive the form of the soft gluon interactions for this situation we can construct an 
auxiliary intermediate theory where the ucollinear gluons and heavy quarks are further from 
their mass shell and the soft interactions are local. The form of this theory is identical 
to Eqs. ( |33||36l , |38"D , but with Y in Eq. fl34|) replaced by QA/m, and we can make the field 
redefinition of Eq. (|39]) in this theory. Then we lower the offshellness of the ucollinear 
particles and match onto the bHQET with scaling exactly as in Eq. (|34|). (This is identical 
to the procedure used to construct SCETn operators from SCETi which was devised in 
Refs. E2j.) The result of this procedure is exactly Eqs. ( |3"BD and (flip. Thus, the result 
for A ^> mT/Q is the same as for A ~ mT/Q. 

We conclude that at leading order the interaction of the bHQET heavy quarks with soft 
gluons are described by Eq. (|4l|) . This matrix element can be used to define a soft function S, 
that describes the cross-talk between massive top-quarks which have fluctuations below the 
mass scale m, and we can use Eq. ([33]) for the remaining dynamics at LO. Thus, the dynamics 
separates in the manner shown in Fig. §, into two decoupled HQET's and a decoupled soft- 
sector. In Sec. [HI F| below we will derive the same result in an alternative manner, starting 
from the factorization theorem for the cross-section in SCET. In this approach the definition 
of the jet functions and soft-cross talk matrix elements are first defined in SCET, and then 
matched onto bHQET. In this case the soft couplings are fully formulated by the matrix 
element in Eq. ([II]), and there is no need to consider soft couplings to fields in the bHQET 
Lagrangian. 

Decay Product Interactions. It is conspicuous that in the leading order bHQET setup, 
gluon exchange involving top and antitop decay products is not present. We now show that 
this treatment is correct and discuss the size of possible power corrections. Since we are 
interested in top/antitop invariant masses in the peak region at large Q, we only have to 
consider ucollinear and soft gluons. Concerning ucollinear gluons it is convenient to switch 
for each bHQET into the respective heavy quark rest frame where v± = (1,0,0,0) and 
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k ~ T « m 

FIG. 4: Example of the cancellation of soft gluon attachments to the decay products. 

the ucollinear gluons have momenta ~ T <C m. For the hemisphere invariant masses 
we can treat the top decay as fully inclusive at leading order (see Sec. PH) , so we can 
address the issue by analyzing possible cuts from the top/antitop final states in electroweak 
diagrams contributing to the bHQET matching conditions [HJ. At leading order in the 
expansion in Y/m there are cuts from the top/antitop self energy which lead to the width 
terms in C±. Subleading finite lifetime corrections to the heavy quark bilinear terms are 
suppressed by Y/m and physically related to the lifetime-dilations coming from residual 
momentum fluctuations of the heavy quark. Furthermore, due to gauge invariance finite 
lifetime matching contributions can not arise for the v± ■ A± couplings in the covariant 
derivatives of £±. Diagrammatically this involves a cancellation between the graphs in 
Fig. |4| including all possible cuts. Diagram a) is a vertex correction, while diagram b), c) 
is a wave-function-type contribution. Since momenta in the cut graphs are of order m, at 
leading order we can take the ucollinear gluons to have momentum = 0. In this situation 
the diagrams cancel due to gauge invariance. Thus, at leading order there are no finite 
lifetime effects involving ucollinear gluon exchange. Effects from the sum of the diagrams in 
Fig. ^ that do not cancel are suppressed by at least a factor a s Y/m relative to the leading 
order factorization theorem. 

Finally we consider soft gluon interactions. Using the proof above for the universality 
of the soft cross-talk matrix element in Eq. fl4"I|) and repeating the arguments made for the 
ucollinear gluon interactions we find that the dominant soft gluon interactions involving 
top/antitop decay products are described by possible cuts of electroweak matching contri- 
butions of the n ■ A s and n ■ A s couplings in Eq. In this case the same cancellation as 
for the ucollinear gluons takes place since the average soft gluon energy in the top/antitop 
rest frame is still A and thus much smaller than m. Thus interactions involving top/antitop 
decays products and soft gluons are suppressed by at least a factor A/m. Numerical studies 
in Ref. @ have estimated QCD interconnection effects based on nonperturbative models 



(see also Ref. |p4| ) 



Having defined the EFT's we now turn to the derivation of the factorization theorem. 
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III. FACTORIZED CROSS-SECTION AND INVARIANT MASS DEFINITIONS 



A. The QCD Cross-Section 



e + e 



We start with the general expression of the cross-section for top-antitop quark production, 
■> 7*, Z* — > ti + X. The final state we are interested in is observed as the top and 
antitop jets plus soft radiation J(t)J(t)X s . We remind the reader that we refer to all 
the jets coming from the top and antitop quark decay collectively as top and antitop jets, 
respectively But we stress that despite the language, our analysis is still perfectly consistent 
with the fact that the different jets from each the top and antitop decay can be resolved in 
the experimental analysis. 
The full cross-section is 



a 



J2WS\q- Px ) £ L l (0|^ f (0)|X)(X|^(0)|0) 



(42) 



x 



where the initial state total leptonic momentum is q = p e - +p e +, Q 2 = q 2 , and the QCD 
currents J^ a are given in Eqs. (|16|). The superscript res. on the summation symbol denotes 
a restriction on the sum over final states X, to give J(t)J(t)X s . These final states contain 
top and antitop jets with invariant masses close to the top quark mass. The explicit form 
of these restrictions depends on the specific jet and invariant mass definitions used. For the 
hemisphere invariant mass prescription these restrictions will be implemented explicitly in 
Sec. piE| below, while other methods are discussed in Sec. [V|. 

In Eq. ( fl2|) we include photon and Z boson exchange, and imply an angular average of the 
leptonic tensor, to obtain the parity conserving U with a sum over vector and axial- vector 
parts, % — v, a. For convenience we also include the charges and boson propagators, and the 
cross-section prefactor 1/(2Q 2 ), so that 
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Here e t is the top- quark charge, and 
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2 sin 6 W cos 9 
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2 sin 9 W cos 9w 



(43) 



(44) 



where T/ is the third component of weak isospin, and 9w is the weak mixing angle. 



B. The SCET Cross-Section 



We now proceed by using the fact that the states are restricted to be dijet-like through 
the constraint that the top and antitop jet invariant masses are close to the top quark 
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mass, as illustrated in Fig. |3|. In this section we reformulate the cross section by using the 
more specific SCET currents of Eq. fll8|) that are suitable for this kinematic situation. We 
integrate out the hard production energy scale Q by matching the SCET currents onto the 
QCD currents giving us via the matching relation ( |T7D a new expression for the cross-section 
defined with matrix elements in SCET. 

The SCET currents in Eq. flIB|) correctly reproduce the long distance physics of the 
QCD current, and the difference in the short distance physics is contained in the Wilson 
coefficient C(u,u, fi). We will see momentarily that momentum conservation dictates that 



the final form of the cross-section depends only on C(Q, —Q, fi) = C(Q, //). In Ref. |65| the 
Wilson coefficient at one loop was computed. It is independent of the Dirac structure Tj and 
also of whether or not the collinear quarks are massive (the latter fact is demonstrated in 
Ref. |)6| where the matching computation for the corresponding vertex diagrams is carried 



out explicitly for finite heavy quark mass). The result is 



OL s Gt, 



-Q 2 -iO , 2 -Q 2 -i0 o 7i 



C(Q,fjL) = l + ^-^ 31og^ V ^4 8 + — . (45) 

At the matching scale [i = Q this Wilson coefficient does not contain any large logarithms. 
The product of the Wilson coefficient C(Q, //) and the SCET matrix element is independent 
of the scale /x, and renormalization group (RG) evolution determines the Wilson coefficient 
at a lower scale /i. This RG evolution of the hard Wilson coefficient sums logarithms of 
fi/Q with /a > m. The Wilson coefficient contains an imaginary part that arises from real 
QCD intermediates states in the QCD vertex diagram that are not accounted for in the 
corresponding SCET diagrams when the collinear action only contains the two sectors for 
the n and n directions (see Sec. |II Aj) . However, only \C(Q,fi)\ 2 will appear in the final 
factorization theorem since we will sum over n. 

Using Eqs. (p~7|) and ( ^8|) in Eq . fl42|) , the cross-section in SCET takes the form 

res. „ 

^ = E E {2-nf5\q-P Xn -P Xf -P Xs )Y J Lf u jdwdwdw'dw' (46) 

n X n XfiX s i 

xC(u,,^)<7Vy)(0|T{x^ 

Here we have pulled out the explicit sum over the top jet label directions n and keep only 
two collinear sectors and for the SCET description of top and antitop jets. This 



allows us to explicitly carry out the integral over the top jet directions n in Sec. [HID 
parallel to implementing factorization. 



m 



In Eq. (f46"D we have decomposed the final states \X) into a soft sector \X S ) and collinear 
sectors \X n ), \X n ) in the n and n directions respectively 

\X) = \X n X n X s ) . (47) 

Since the hard production scale is integrated out by the matching procedure, these states 
now form a complete set of final states that can be produced by the SCET currents . 
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This already implements part of the restrictions, "res", in the sum over states in Eq. (fig). 
The momentum Px of the final state \X) is also decomposed into the momentum of the 
collinear and soft sectors: 



x 



P Xn + Px n + Px s . (48) 



Recall that there are no particles with pf^ ~ (m, m, m) scaling that can cross the final state 
cut, without taking the invariant mass far from the peak region, so there are no mass- modes 
in this decomposition. Because the set of hadrons observed in the detector has a well defined 
set of momenta, it is possible to impose criteria on the hadrons in the final state to associate 
them with one of X n , X n , or X s . Thus, the hadronic two-jet state factorizes as a direct 
product 

\X) = \X n )\X n )\X s ) . (49) 

This factorization is also a manifest property of the hadronic states in SCET. 

For quark and gluon states in SCET the difference from the purely hadronic case in 
Eq. ( [4"9"D is that the states can carry global color quantum numbers. After having made 
the soft-collinear decoupling field redefinition, the individual Lagrangians for these sectors 
are decoupled, and they only organize themselves into color singlets in the matrix elements 
which appear in the observable cross-section. We can take this as a manifestation of quark- 
hadron duality. Using the soft-collinear decoupling property from Sec. [11 A| we can write the 
matrix elements in Eq. (HQ) as M. (m, fx) times 



(ofe^) 6a ^(ry nW Ol^nX s x s ><x n x s x s |(x n ,^r) c (yt)^|o) (50) 

= <o|x^|x,)<x,|x^|o><o|x^|x n )<x n |xlJo> 
x (o|(F n ) co (ry n ) cfe |x s >(x s |(y n tr) 6 ' c '(Fi) £lV |o) , 

where here roman indices are for color and spin and \X n ) and \X n ) are color triplets. Next 
we rearrange the color and spinor indices so that they are fully contracted within each of the 
n-collinear, n-collinear, and soft product of matrix elements. This makes explicit the fact 
that in SCET each of these contributions to the cross-section must separately be a spin and 
color singlet. Although it is not absolutely necessary to make this arrangement of indices 
manifest at this point, it does allow us to avoid carrying around unnecessary indices (a similar 
manipulation was used for B — > X s ^ in Ref. j67|| ). For color, our \X n ){X n \ forces the indices 
on xl and xt to be the same, so (o|^|X ft )(X ft |x£|0) = (S aa ' /N c ) (o|*> |X fi >(X fi |x* |0). A 
similar result holds for the n-collinear matrix elements. For spin we can use the SCET Fierz 
formula 



1 ® 1 



1 



2 L v 2 



(f)^) + rf>#)+(4%#)l. (si) 
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which is valid when the identity matrices are inserted so that the t/i terms on the RHS 
appear between Xn" 'Xn without additional fi factors next to these fields (or the analogous 
statement with n n). Combining the color and spin index rearrangement, the matrix 
element in Eq. (BUT) becomes 



tr 



L2 i 2 j . 



°>1 \(°\(i^rXny) b \x n )(x n \t n Jo) 



x 



4N c Xn 'V '"'J L^'V4iV c 



= tr 



<o|(r,) ca '(F n ) cfe '|x s ><x s |(F n t) 6V (rI) a ' c '|o> 

tr(<0|x^|^)<X s |^Xn,cD|0>)tr(<0|^Xn^|X n )<X n |x n , w |0> 



xtr^(o|y fi y n |x 8 )(x a |iJyi|o) 

where for convenience we defined 



ft = i/{±N c 



fL = i/{AN c ) . 



(52) 



(53) 



Note that only the first term on the RHS of Eq. ([H]) contributes because the collinear states 
give at least one matrix element which is zero when we have a 75 or 7". This factorizes the 
SCET cross-section into a product of three singlets under spin and color. For convenience 
we will in the following suppress writing these explicit traces on the matrix elements. 
Using Eq. fl52j) in Eq. fl46|) , the factorized SCET cross section takes the form 



a =K MJ2 E (^) 4 S 4 (q-P Xn -Px,-Px a )(0\Y n Y n \X s )(X s \Y^Y{\0) 

71 X n Xfi X s 



(54) 



Xj duo duo duo' duo' C(iU, u)& (u', uo') (0\jtXnV \X n ) (Xn\Xn,JO) {^XnM X n) { X nWn,M 

where M. = M. (m, fi) and we defined the normalization factor 



V L«Tr 



2 1 2 



32n 2 a 2 



= -2gT E L S 

i=v,a 

2Q 2 v e v t e t , Q 4 (v 2 + a 2 e )(v 2 + a 2 ) 



Q 2 — m| 



+ 



(Q 2 - m|) 2 



(55) 



We can further simplify the form of the factorized cross-section. First we use the identities 



{ X n\Xn,M = ( X n|Xn^',n-Pt|0) = ^, ip -^ (X n |x n |0), 
( X n\Xn,oji\ty = ( X n\xju',n-V^) = 6 -u',p+_ ( X n\Xn\fy 



(56) 

with similar relations for the other two collinear matrix elements in Eq.(j54"D. Combining this 
with the relation 8 W , - 5^ - = - , and analog for p~^_ , we can write the product of 

collinear matrix elements in Eq . (|54]) as 

(0|fenxl^n)(^n|X^|0)(0|Xn^l^n)(^n|fea)|0> 

= 6^6^ (0|fe|^n)(^n|X„, w |0)(0|Xnl^n)(^n|fea)|0) • (57) 
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Next we do the sums over u', u' to arrive at the form 

res. 

a = K MJ2 E ^)U\q-P Xn -P Xn -P Xs ^\YnY n \X s )(X s \Y^Yl\Q) 

n X n Xn X s 

xJdu;du\C(uj,oo)\^0\fiXn\X n )(X n \x n JO)(0\xn\X^ (58) 

Before proceeding, we pause to define the thrust axis which is needed to properly define 
the invariant mass of jets and to state its relation to the direction of the energetic collinear 
degrees of freedom. Then in order to make the power counting manifest we decompose the 
final state momenta into label and residual parts and perform some general manipulations of 
the phase space integrals to setup a formula for the cross-section to be used for the remaining 
calculation. 



C. Thrust or Jet Axis 



The thrust T of any event is defined to be 

T = max ^ , (59) 

t Q 

where the sum is over the momenta of all the final state particles produced. The thrust 
axis t is chosen so that is maximizes the sum of particle momenta projected along t. Intu- 
itively, for a dijet-like event the thrust axis corresponds to the axis along which most of the 
momentum is deposited. Conversely, the thrust T is close to its maximum for a dijet-like 
event. We choose n to point along t. For an event with exactly two massive stable particles 
T = ^Q 2 -Am 2 /Q = 1 - 2m 2 /Q 2 + 0(m 4 /Q 4 ), is the maximum allowed thrust. Since we 
are interested in thrusts in the dijet region for the top and antitop jets it is convenient to 
define a shifted thrust parameter, 



4m 2 _ 2m 2 m M / m 4 \ 

1 -^- r = 1 -^- r+ °y- < 60 > 

For stable top-antitop production additional jets always result in r > 0. For unstable top- 
quarks the values of r < also become allowed. Note that for massless jet production the 
thrust (T) distribution is peaked close to T = 1 while for events containing a heavy quark 
pair it is peaked close to T = \/Q 2 — Am 2 /Q. Thus a cut on thrust can in principle be used 
to discriminate between massive and massless quark production ||. 

D. Differential Cross-Section with Momentum Decomposition 

To insert the invariant mass constraints into our cross-section in Eq. ([54]) we use the 
identity operator: 

1 = / d 4 Pn d% d A Vs 5\ Pn - P Xn ) 6\ Pn - P Xn ) 6\ Ps -P Xs ), (61) 
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which sets the total collinear and soft momenta of the states Px n , Px n , Px s to p n} Pn,Ps 
respectively. In Sec. E| we will use an additional insertion of an identity operator to define 



the hemisphere invariant masses, M t and Mj. In this section we carry out manipulations 
that are common to any definition of the invariant masses. For now we ensure that the 
invariant mass of each hemisphere is close to the top mass by including in the restrictions, 
"res", on the states the fact that M t , Mj are in the region 

I >»*,*| = \ M t.t - m 2 | < m 2 . (62) 

From here on we assume that in the sense of power counting A ~ V. We now decompose 
the collinear and soft momenta into label and residual parts 

Pn = Pn + k n , Pn = Pn ~\~ kft, ^J(„ = K x n i (63) 

^Xfi ^ X n ' ^ X n Xn ^ X n ' ~^X n ^X n ^ X n ' 

^X n ^ X n ' ^ Xfi ^ ~ Xn ' ^*X S ^ X s ' Ps kg 

Note that our choice of n along the thrust axis together with the restrictions on the states 
ensures that the perpendicular momentum of the jets relative to this axis, P% and P^, are 
purely residual. The last result in Eq. fl6"3| ) indicates that the soft state also has a momentum 
that is purely residual. The integrals in Eq.([H]) can be decomposed into a sum over labels 
and integrals over residual momenta as 

J d'p n J d*p n = \Ylj d KdKd 2 ki \Ylj dk t dk nd 2 ki . (64) 

pn Pn 

In the total cross-section in Eq. fl54|) we sum over the directions n of the thrust axis. To 
turn this sum into an integral over the full solid angle, we need to combine it with a residual 
solid angle integration for each n. Therefore, we decompose the residual measure as 

(Q 2 \ 

d 2 kn — \p n \ 2 d(j) d cos(9 r ) = I— p 2 n )d(f) dcos(9 r ), (65) 



4 

where 8 r is the small angle of p^ relative to p. In the first equality we used the fact that 
cos(6* r ) ~ 1. Since we are in the peak region we can approximate p 2 n = m 2 up to small T/m 
corrections. Combining this with the sum over n gives 

rG 2 \ O 2 
J^d 2 k^ = (-j- -m 2 ) d0dcos(9) = -^-dtt, (66) 

n 

where in the last equality we work to leading order in m 2 /Q 2 . Since the angular averaged two- 
jet production is independent of the thrust direction we are free to carry out the remaining 
integrations in a frame where = 0, and also replace J dQ = An. The differential cross- 
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section now reads 

a= ^_J> M J2 ( 2 n) 4 5 4 (q-P Xn -P Xn -P Xa ) J2 JdktdK J dk+ dk~ n d 2 % d% 

X n XfiX s pn-iPn 

x 5\ Pn - P Xn ) 5\p n - P Xn ) 5\k s - P Xs ) (Q\Y n Y n \X s )(X a \Y^Yl\U) 

x ^0;^ |C'(^,^)| 2 <0|^(0)|^><^|x^(0)|0><0|xfi(0)|^><^|fe^(0)|0> . (67) 

In the remainder of this section we will simplify this formula as much as possible prior to 
specifying the exact constraints on the restricted sum of states. First we decompose the 
delta functions into label and residual parts as 

£4/ _p s _ 5 . 5 4 (k _ R ) _ 6 g f_^_ i[(ki-K+J^ + (^-K Xn )^-K^. X ^ 
\Pn ^X n ) - °p n ,P Xn l ft n J^X n ) - J J > 

t\Pn-PxJ = 5^ Px j\k n -K x J = 5^_J Ph0 J-0y -v, 

5\Vs-Px s ) = 5\k s -K Xs ) = Jj0ji e^~ K ^-\ (68) 

where there is no k^ in the first line (or below) because we fixed k^ = 0. In the second 
equality on lines 1 and 2 we replaced P^ n ,P^_ with the labels u,u respectively using the 
momentum conservation delta-functions discussed below Eq. (|56|). We also decompose 



5\q-P Xn -P Xn -K Xs ) = 5 QS - 5 QtPt 5\k n + k n +k s ) , (69) 

where we have replaced P Xn , P Xn with p n ,Pn by using the delta functions in Eq. (|68|). 

Next we use Eqs. ( j68|) and ( |69"D in Eq. (p7| ) and with the exponential factors of 
e~ tKx "' x ,e~' lKx "' y , and e~ lKxs ' z in Eq. ( p8[ ) we translate the collinear and soft fields to the 
positions x, y, and z, respectively. This gives 



vr Q 2 K 

a 



(2vr) 8 4 



Ai Jdk^dk n dk^dk n d 2 k^d 4 k s Jd 4 x d 4 y d 4 z 5 4 (k n + kn + k s 



X n X n X s 



x \C(Q,fi)\ 2 ExpJ-fc+aT + l -k~x + + ik n -y+ik s -z 1 \ (o\(Y n Y n )(z)\X s )(X s \(Y r lY j n )(0)\0) 
x (0|^Xn(^)|^n><^|x„,o(0)|0>(0|x ft (y)|X fi >(X fi |te.-o(0)|0) , (70) 

where here the large label momenta in the jets are fixed to be Q, X n ,Q = Xn^Q,v^ an d 
Xn-Q = 8-Q,vXn- Next we can use the fact that the n-collinear graphs are independent of the 
k~ and k^, so that the above n-collinear matrix element is proportional to S(x + )S(x±) |28| . 
Similarly the n-collinear matrix element is oc 5(y~)5(y±). It is not crucial to use these 5- 
functions at this stage, but they do allow us to simplify the formula by dropping x + , x±, 
y~, and y± dependence in the exponentials. Performing a few integrals we arrive at a fairly 
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simple form for the cross-section 

a = a \C(Q,fi)\ 2 M Jdk+ dfcr dk+ dk; (71) 

res. 1 „ 

Jd 4 xe ik ^-/ 2 (0\fiXn(x)\X n )(X n \x n>Q (0)\0) 



res 

X 

X 

res 

X 



JdSe^y + ' 2 {Q\xM\Xn){Xn\fiXn,-Q{^) 

res. 1 „ 

x Ei]v^ Jdz+dz-e^*- +k °^ (0\Y n Y n (z-,z+)\X s )(X s \YiYl(0)\0) . 
x s 

The result in Eq. (|71|) is a factorized product of Fourier transforms over n-collinear, n- 
collinear, and soft matrix elements. We introduced a 1/N C in front of the soft-matrix element 
in Eq. ([TJ]), and include a compensating factor iV c in cr - This equation provides a good 
starting point for the derivation of any differential cross-section (for massive or massless 
dijet events). The new normalization factor <7 is just the total Born cross-section 



Q 2 „ Ar 4ttc^ 



vr u c 3Q 



2Q 2 v e v t e t , Q 4 (t; e 2 + a2)K 2 + a?) 



g 2 -m| (g 2 -m 2 ^ 2 



(72) 



For massive quarks <To depends on (3 m = (1 — Am 2 / Q 2 ) 1 ^ 2 through an extra multiplicative 
factor of /3 m (3 — /3^)/2 = 1 — 6m 4 /Q 4 + . . .. This is only a 1% correction to a for (J/m ~ 5. 

To proceed further we now need to make the prescription how the n- and ra-collinear and 
the soft particles enter the invariant masses St and sj explicit. This removes the implicit 
restrictions on the sums over states indicated in Eq. ([71]). In the next subsection we imple- 
ment the prescriptions for the hemisphere jet invariant masses. In Sec. [V] we briefly discuss 
how the implementation changes for other invariant mass prescriptions. 



E. Factorization for Hemisphere invariant masses in SCET 



In the hemisphere mass case all the final state particles are assigned to belong to one of 
two hemispheres defined with respect to the thrust axis. The boundary between the two 
hemispheres is perpendicular to the thrust axis and centered at the e + e~ collision point, 
see Fig. |^. Thus the top and antitop jets we consider correspond to all the particles in 
the respective two hemispheres and the invariant mass of each jet is defined to be the total 
invariant mass of all the final state particles in each hemisphere. As we show explicitly below, 
the requirement that these jet invariant masses are both close to the top mass, automatically 
restricts the final state to be dijet-like, and eliminates the need to introduce any additional 
event-shape constraint. We stress that some mechanism to control the soft particles is 
absolutely crucial for establishing the factorization theorem and the unique definition of 
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the soft function S. Here this is accomplished by the fact that all soft particles enter the 
invariant mass variables M 2 ^. 

The invariant mass of each hemisphere includes contributions from both soft and collinear 
particles. The total momentum of the collinear particles in the n-hemisphere is Px n and in 
the n-hemisphere is Px n - The total final state soft momentum Kx s is split between the two 
hemispheres and can be divided as: 

K Xs = K + k) (73) 

where and k b correspond to the total momenta of all the soft partons in the n and n 
hemispheres, respectively. It is useful to think of these hemisphere momenta as the result 
of hemisphere projection operators P a , Pt,: 

P a \X s ) = k a s \X s ), P b \X S ) = k b s \X„). (74) 

In other words, these projection operators act on each state \X S ), pick out the soft partons in 
the respective hemisphere and add up their total momentum. Note that the eigenvalues are 
dependent on the state X s , so k" = kg[X s ] and k b s = k b [X s }. We can now define the invariant 
mass of each jet as (Px n + K) 2 and (Px n + k b s ) 2 for the n and n hemispheres respectively. 
The delta functions 5 A {p n — Px n ) <5 4 G°n — Px n ) m the second line of Eq. (|67D allow us to 
define the jet invariant masses in terms of p n ,Pn as (p n + k^) 2 and (p n + k b s ) 2 for the n and 
n hemispheres respectively. 

Note that this implements a very simple form of a jet algorithm. For a different jet 
algorithm we would change the definitions of the operators P a and Running a jet 



algorithm in inclusive e + e~ mode |68, each soft parton is still accounted for, having a 
certain probability of being assigned to either the top or the antitop invariant mass. We 
discuss other algorithms in Sec. 0. 

If the top quark were a stable particle, these invariant mass definitions would be obvious 
because n- and n-collinear particles would be fully radiated into the n- and n-hemispheres, 
respectively. Due to the finite lifetime of the top quark, however, we need to convince 
ourselves that this invariant mass definition still works if the n- and ra-collinear momenta 
of the top and antitop quarks, respectively, are distributed among their decay products. So 
let us consider the top quark in the n-hemisphere. Since the top rest frame is boosted with 
respect to the e + e~ cm. frame with a boost factor Q/m, top decay events can have final 
state particles appearing in the n-hemisphere of the antitop quark only if these final state 
particles have an angle (defined in the top rest frame) smaller than m/Q with respect to 
the antiparticle direction. On the other hand, the top spin is only about 20% polarized 



(for unpolarized e + e~ beams and upon averaging over the directions of the thrust axis) [70 
and thus the top decay products in the top rest frame are distributed isotropically to a 
rather good approximation. The fraction of events in this kinematical situation is therefore 
suppressed by (m/Q) 2 and can be neglected at leading order in the power counting. Of 



30 



course the analogous conclusion also applies to the antitop quark in the n hemisphere. So at 
leading order in the power counting it is consistent to employ the invariant mass definition 
of the previous paragraph. 

The jet invariant mass definitions can be implemented into the cross-section of Eq. ( |7T| ) 
by inserting underneath the Ylx ^he identity relation 



dM 2 t 5((p n + Kf - Ml) J dMl 5(( Pn + k b ) 2 - Ml) 
= JdM? 5((p n + k a s ) 2 - m 2 - s t ) J dMl 5((p n + k b s f - m 2 - s f ), (75) 

where s t = s t (M t ) and sj = sj(Mf) from Eq. ([!]), i.e. it should be understood that s t j are 
functions of M^. In the second line m is defined as the pole mass. It is straightforward to 
switch the final result to a suitable short distance mass definition, as we explain in Sec. [Ill G . 
Decomposing the 5-functions at leading order gives 

8((p n + Kf -m 2 -s t ) = ± 5{k + n + kf a - ^±^) , 
6((p n + k b ) 2 -m 2 - at) = - 5(k7 n + k 7 b - !^±5) , (76) 



Q \ n 3 Q 

where we set p~ = p£ = Q due to (^-functions from Eq. (|69|) . Carrying out the integration 
over fc+ and k~ in Eq. ( |7T| ) sets the arguments of the soft function to z* 1 = 0. Inserting the 
identity relation 

'd£ + d£-5(£ + - kj a )5(£- - k~ b ) (77) 



the differential cross-section then reads 



d ^f dM , = ^\C(Q,fi) 1 2 M J ikt dk~ n d£ + d£~5(k+ +£+- 5 (k~ n + 
£ A. ke'^' 2 tr(0\fiXn(x)\X n )(X n \x n , Q (0)\0) 

V ^ J 



t 

x 



X 



JdSe^/ 2 tT(0\xn(y)\Xn)(X n \fi X n,-Q(0)\0) 
J2w 6 ( £+ ~ ~ k7 b )tr(0\Y n Y n (0)\X s )(X s \Y^YU0)\0) , (78) 



X 



where we have dropped the "res." label on the sums, because all restrictions are now 
explicitly implemented. 

To see that the hemisphere definition of s t and sj can be used to select dijet-like events, 
we can check that Eq. ( j62f) plus the 5-functions in Eq. ( j78|) are sufficient to constrain the 
thrust to the dijet region. At leading order the total thrust of an event is given by 

Q T = bnl + \Pn\ + + \ k a*\ j (79) 
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where 2|p*| = Q + k n — k£ and 2|p?| = Q + k£ — k n . So the shifted thrust defined in Eq. (|60] 
is 



2m 2 

Q 2 + 


1 

2Q 


2m 2 

Q 2 + 


1 " 
Q- 


St + St 




Q 2 ' 





h + — h~ -I- h~ — k + 4- h a+ — h a ~ 4- k°~ — h 

k^ -\- k -\- -\- k„ 



(80) 



To obtain the second line we used the separate conservation of the + and — momentum 
components to eliminate k~ + k b s ~ and k% + k^ + . For the last line we used the ^-functions 
in Eq. ( [TED to get s t and si- Thus, the restriction to small hemisphere invariant masses s ty i 
automatically gives small r and restricts the events to the dijet region. The presence of a 
third hard jet takes us away from the dijet region and directly shows up by a substantial 
positive shift of s t + away from the peak region. 

Next we simplify the form of the cross-section by defining the massive collinear jet func- 
tions J n , J n as 



/j4 
^ e -^J n (Qr+-m 2 ,m) (81) 

= Q5(x + )5 2 (x ± ) Jdr+ e-* ri *- J n {Qr+ - m 2 , m) , 
J2te(U\Xn(v)\Xn)(Xn\fiXn,-Q(0)\0) = Q J ^e~^ y J n (Qr^ - m 2 , m) 

= Q5(y-)5 2 (y ± ) J dr^ e ~^ y+ J n {Qr^ - m 2 , m) . 

Here m is the pole mass just as in Eq. (|75|) and we do not display the jjt dependence. 
Note that the subscript Q on the LHS does not change the mass-dimension of a x-field 
away from 3/2, since 5np is dimensionless. We remind the reader that fi = rf/(4N c ), p 
= t/i/ (4iV c ) and that tr is a trace over both color and spin indices. The arguments of the 
jet functions, J n and J n , in Eq. fl8~T|) are just the off-shellness of the jets, p 2 — m 2 and 
— m 2 , respectively, but given in expanded form. Here the labels Q on the x n an d Xn 
fields ensures that there is only a contribution from the required "quark" and "antiquark" 
cut since Q > 0. To see this recall that the sign of the label p on picks out the 
quark annihilation, or antiquark production part of the field ||27|| . We note that the sums 



over collinear states in the collinear jet functions are unrestricted since the restrictions are 
now implemented automatically through the amount the jet invariant mass differs from m 2 . 
Thus, the jet functions can be written as the discontinuity of a forward scattering amplitude 
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after summing over the collinear states: 

J n {Qr+-m\m) = -±Discjd 4 x e ir - x (OiT{x n , Q (0)fc(^)}|0) , 

J- n {Qr- n -m\m) = -L Disc [d 4 x (0\T{xn(x)fi X n,- Q (0)}\0) • (82) 

ZttQ J 

The collinear fields in the SCET jet functions J n and J n are defined with zero-bin sub- 
tractions [54], which avoids double counting with the soft-function. Using Eq.(|8T|) and 
performing all the remaining integrals in the cross-section of Eq.(^) we arrive at the SCET 
result for double differential hemisphere invariant mass cross-section 

a H Q (Q,fi)M(m,fi) (83) 



x / d£ + d£ J n (st-Q£ + ,m,ij,)J n (st-Q£ ,m, fi)S hemi (£ + ,£ ,/i,m), 

J — oo 

where the hard function Hq(Q, fi) = \C(Q, /u)| 2 . Here the hemisphere soft function in SCET 
is defined by 

S hem i(^,r,//,m) = -L Y, 5 ^ + - k t a W ~ K b )(0\YnY n (0)\X s }(X s \Y^Yi(0)\0) . (84) 

At tree level for stable top quarks Hq = 1, J n (s t ) = S(s t ), Jn(st) = S(st), and S , h em i(^ + , £~) = 
S(£ + )5(£~), and integrating Eq. ( |8"B| ) over s f and Sf gives the total tree-level Born cross-section 
o"o- This provides a check for the normalization of Eq. (|8"3]). The argument m on the soft- 
function in Eq. (0) and A4(m, /j,) in Eq. (|83D account for massive top-quark bubbles that 
are perturbative and start at 0(a 2 s {m)) |7l], [72], |73|. Note that Eq. (^) extends the SCET 



computation of the massless dijet cross-section in Ref. |34], |57| to all orders in perturbation 
theory for the jet-functions. 

In the factorization theorem in Eq. (|83|) the jet-functions J n and J n describe the dynamics 
of the top and antitop jets. In the next section we will see that these jet functions can be 
computed in perturbation theory and at the tree level are just Breit-Wigner distributions. 
The soft matrix elements (0\Y^Y n (0)\X s ) (X s \Y^Y n (0)\0) , on the other hand, depends on the 
scale Aqcd, and thus the soft function Sh em i(^ + , £~) is governed by non-perturbative QCD 
effects. The momentum variables £ represent the light cone momentum of the soft particles 
in each of the two hemispheres, and Sh em i(^ + , £~) describes the distribution of soft final state 
radiation. Eq. ( p3|) already demonstrates that the invariant mass spectrum for unstable top 
quarks is not a Breit-Wigner function even at tree level because the convolution with the 
soft function Sh em i modifies the observed distribution. The effects of the convolution on the 
observable invariant mass distribution are discussed in Sec. [TV]. 

To sum large logs in Eq. (|83|) the SCET production current can be run from /j, = Q 
down to fi = m, which then characterizes the typical virtuality of the collinear degrees of 
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freedom in massive SCET. In the process, large logarithms of Q/m are summed into the 
hard function Hq(Q,(i). In the next section we integrate out the scale m and match these 
SCET jet functions onto bHQET jet functions. 



F. Factorization of Jet mass effects in HQET 

The main result of the last subsection is the factorization of the scales Q and m in the 
differential cross section of Eq. fl8~3"|). In this section we further factorize the scale m from the 
low energy scales T, s, and A. This will allow us to sum large logs of T/m and s t j/m in the 
jet functions, and lower the scale of the soft functions to A. This step is also important for 
treating the width effects. As explained earlier, one can formulate width effects in a gauge 
invariant way with a natural power counting in HQET, whereas doing so in a relativistic 
theory such as SCET is notoriously difficult. 

To perform the scale separation and sum the logarithms requires us to match and run 
below the scale /i = m. This can be done in a standard way, by matching and running of the 



bHQET current in Eq. (p7|), as we described in Sec. [II B| . However, due to the factorization 
properties of SCET which leads to a decoupling of the n-collinear, n-collinear, and soft 
sectors, the matching and running below the scale [i = m can also be done independently 
for J n , J n , and S. In the following we explain this second method. 

As discussed in Sec. [IIB| the soft function above and below the scale m is identical, except 
for certain vacuum polarization effects from graphs with top-quark bubbles that only exist 
in SCET. For the soft-function in bHQET we have 

5 hcmi (^,r,/i) = 1- - K a W - K b )(0\Y n Y n (0)\X s )(X s \Y^Yl(0)\0) , (85) 

c x s 

where there is no m dependence. The matching condition between the soft functions in the 
two theories is 

Shcmi(£ + , m) = T (m, /i)S hemi (£ + , £~,fi), (86) 

where To(m, fi) is a Wilson coefficient. Large logarithms in the soft function can be summed 
by computing the anomalous dimension of the soft function and using RG evolution to run 
between A and m and between m and Q as illustrated by the line labeled Us in Fig. |5|. 

For the SCET collinear degrees of freedom the power counting for the virtuality is p 2 ~ 
m 2 . Thus, J n and J n describe the physics of jets with an invariant mass up to M 2 ~ fi 2 ~ m 2 . 
However, the restriction of being in the peak region means that M 2 — m 2 <ti m 2 . This 
disparity gives rise to the large logarithms in the collinear jet functions. Intuitively, this can 
also be understood by noting that if one starts out with a top quark that is close to its mass 
shell, a typical collinear SCET gluon will knock the top far offshell so that p 2 — m 2 ~ m 2 . 
By restricting the jet functions to p 2 — m 2 m 2 we forbid such real radiation contributions, 



34 



Matching Functions 
Hq(Q, M 



H m (m, IV) 
S(A, |I A ) 

B + (r, |i r ) 



Top-down 
running 



Bottom-up 
running 



Scales 



> 


> 

U Ho U Jn 

< 


k > 


\ / 






> 

u Hm 

U B . 

t 


k 1 


k / 
U B+ 


Us 


U B .|1 U b+ 





m 
A 

s/m ~ r 

Aqcd 



FIG. 5: Scales and functions appearing in the formula for the invariant mass distribution. The 
result is determined by matching at the physical scales and running to sum large logs as shown. We 
show both the top-down and bottom-up approach to the running. The evolution for Uh and Uc 
is local, while all other evolution functions involve convolutions. Note that the evolution functions 
obey Uh = Uj_ ® Uj + <g> Us and Uc = Ub_ <8> Ub+ Us where (g> indicates convolutions. 

but not virtual contributions. The latter must be integrated out explicitly by switching 
to the description of the jet functions in the boosted unstable HQET theories discussed in 
Sec. [II B| . In these HQETs the only fluctuations are due to low energy ultracollinear gluons 
that preserve the condition M 2 — m 2 <C m 2 . 

To determine the definitions of the bHQET jet functions we follow the same procedure as 
for the bHQET current in Eq. (|37|) , namely boost the SCET jet function in Eq. (|8^) to the 
heavy quark rest frame, giving ip(x)W(x)W(0)ip(0), then match onto HQET ift(x) — > h v (x). 
We then boost back to the moving frame where v — > v±. The spin structure can also be 
simplified to give 



4iV c Q 

Thus the bHQET jet functions are defined as 

-1 



V -±^h h 



4N r m hv+hv+ 



(87) 



B + (2v+-k) 



lnN c m 
1 



Disc jd A xe ikx (0\T{h v+ {0)W n {0)Wl{x)h v+ {x)}\0) 
Disc jd A xe ikx (0\T{h v „(x)W n (x)W^0)h v _(0)}\0). 



(88) 



8nN c m 

These bHQET jet functions can be calculated using the usual Feynman rules of HQET 
except that the gluons have ucollinear scaling as in Eq. (0). The W- Wilson lines in B± 



also contain these boosted gluons. Since p 2 



m 



2mv + ■ k and p 2 



m 



can identify the arguments of the bHQET jet functions as 

St o ; s * 

— = St , 2v- ■ k = — 

m m 



2v+ ■ k 



si- 



2mv- • k, we 



(89) 
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a) 

® > 



b) 



FIG. 6: Tree level top-quark jet functions in a) SCET and b) bHQET. 



In the factorization theorem these arguments are shifted by the soft gluon momenta as shown 
in Eq. (|94]) below. Recall that the fields h v+ and h v _ are defined with zero-bin subtractions 
on their ultracollinear momenta. For Eq. fl88|) these subtractions can be thought of as being 
inherited from the SCET fields in the matching. They remove the light-cone singularities 
n ■ k — > and n ■ k — > in B + and B_ respectively, and are important to ensure that the 
width T is sufficient to make B± infrared finite. 

In general the matching of the jet functions in SCET onto those in bHQET could take 
the form 



where the convolution takes into account the fact that depending on the definition, the 
observable s could be sensitive to scales of 0{m) and 0(T). In such a case, since s' does not 
know about the scale m, it can not be identical to s. The convolution with T±(s, s' , m, fi) 
then compensates for this difference. In our case (and most reasonable cases) the definition 
of the invariant mass is not sensitive to m, so we have T±(s, s' , m, /x) = 8(s — s')T±(m, /i) 
and the matching equations are simply 



Since there are no mass-modes in bHQET the function Ad also appears as part of the Wilson 
coefficient. From this we define a hard-coefficient that contains the mass corrections as 4 



By charge conjugation we know that the jet functions for the top and antitop have the 
same functional form, and that T + = T_. When we sum large logs into the coefficient H m it 
develops an additional dependence on Q/m through its anomalous dimension which depends 
on v+ ■ n = v _ • n = Q/m. Note that in principle H m (m, /i) and the factors in Eq. ( P^) can 
also have Q/m dependence at NNLL order. For related discussions see Refs. JF4J [TEfl . 

4 In explicit computations scheme dependence may effect the manner in which the mass-corrections are 
divided up between T±, Tq, and A4, however this dependence will cancel out in the product that gives 




(90) 



J n (ms, m, T, /i m ) = T + (m, /i m ) B + (s, T, /i m ) , 
J n (ms,m,T,fx m ) = T_(m,/i m ) B_(s,T,n m ) . 



(91) 



H m (m,fx m ) = T + (m,/i m )T_(m,/i m )T (m,/i m )A^(m,/i m ) . 



(92) 



H, 



m • 
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Since the functions T± are independent of the top width T, we are free to set T = (i.e. use 
stable top quarks) for the matching calculations at any order in perturbation theory. At 
tree level we need to compute the discontinuity of the graphs in Fig. |] which have a trace 
over spin and color indices. For T = this gives 

Bf c (s,T = 0) = ~l (—2N C ) Discf \ -) = — *— Im 



wN c m \v + -k + iOJ 47rm \v + ■ k + iO 

= — 5{2v + ■ k) = — 5{s) = 5{s) , (93) 
m m 

which is identical to the result for the corresponding SCET jet function, so at tree level 
T + = T_ = 1. 

Plugging Eq. fl9~l|) into Eq. ([83]), and incorporating renormalization group evolution, the 
form for the differential cross section is 

( dM 2 dM? ) = cr oH Q (Q,fx m )H m (m,^-,fx m ,f?j (94) 

> * t ' hemi 

x / dtdr b + ($ t - r, fx) b_ (h - r, fj) s hcmi (£ + , r,/i). 



Eq. (0) is our final result in terms of the pole mass m. The analogous result for a short 
distance mass is given in the next section. Here H m (m, Q/m, /i m , /i) is the hard coefficient 
H m (m, fjL m ) run down from /i m to /i, and we still have HQ(Q,fi m ) = \C(Q, /i m )| 2 , and the 
soft function with Wilson lines evaluated at x — 0, 



(95) 



c x s 



For completeness we wrote out the color indices from Eq. (|52|). Its interesting to note that in 
the result in Eq. ( [94] ) the final matrix elements only involve Wilson lines (since the coupling 
of gluons to a heavy quark field h v+ in B + is the same as to a Wilson line W v+ ). 

To conclude this section we finally repeat the computation of the tree level bHQET 
jet functions, but now for the realistic case with T 7^ in the HQET propagators. The 
computation is done at a scale > T, but the \i dependence does not show up at tree level. 
Fig. ^|b gives 

Bf e (s,T) = -^—(-2N C ) Discf / — ) = -*-Im( ,~ 2 - , ) 

± v ' ; 8vriV c m v ; \v±-k + iT/2) Anm \v±-k + iT/2J 

1 r 

irm s z + I z 

Thus we see that B±{s) are equal to Breit-Wigners at lowest order in a s . At higher orders 
in perturbation theory the width will cut off the IR divergences that would otherwise occur 
at s — 0. The functions B± at the scale \i can therefore be computed perturbatively 



37 



to any desired order in a s . In general the perturbative "matching" corrections will lead 
to distortions of the tree- level Breit-Wigner distributions shown in Eq. (pop, as does the 
potential separate running between /ia and fir discussed below in Sec. |111H| and shown in 
Fig. § 



G. A Short-Distance Top-Mass for Jets 



The derivation of the factorization formulae (|94|) in the previous section was given in the 
pole mass scheme 5 , m po i e . It is, however, well known that the pole mass definition leads to 



an artificially enhanced sensitivity to small momenta in Feynman diagrams (see Ref. [75| for 
a review) and, as a consequence, to artificially large perturbative corrections. This behavior 
is particularly important for observables that have a strong dependence on the heavy quark 
mass [|1^, |58|, |5§|, |6(], |61| . From a nonperturbative point of view, this feature is related to an 
intrinsic ambiguity in the heavy quark pole mass parameter of order the hadronization scale 
Aqcd, and is sometimes referred to as the 0(AQCD)-renormalon problem of the pole mass. 
Heavy quark mass definitions that do not have such an O(Aqcd) ambiguity are called short- 
distance mass schemes. 6 In the factorization formulae in Eq. (0), the top-mass appears in 
the hard function H m and in the two jet functions B + (s t ) and B_(§t)- The most important 
sensitivity to the top-mass scheme is in s t = [M 2 — m 2 )/m and sj = (Mr — m 2 )/m, where 
M 2 and M 2 are scheme independent observables. 

A specific short-distance top quark mass scheme "m" can be defined by a finite residual 
mass term 5m ^ 0, as 



m 



pole = m + 5m, (97) 



where 5m starts at 0(a s ) or higher, and must be strictly expanded perturbatively to the 
same order as other 0(a s ) corrections. (This strict expansion does not apply to powers of 
a s times logs that are summed up by renormalization group improved perturbation theory.) 
Let B + (s, fi,5m) denote the jet-function in the short-distance mass scheme specified by 
5m. We can calculate B + (s, fi,5m) in two equivalent ways, i) Use the pole- mass scheme 
initially by setting 5m = in Eq. (1331). In this case the mass-dependence appears in 



Spoie = (M 2 — m 2 olc ) I 'm po ic in B + and we change the scheme with Eq. fl9~7l) . Alternatively, 
ii) treat 5m ^ in Eq. fl33|) as a vertex in Feynman diagrams, and take s to be defined in 
the short- distance mass scheme right from the start, so s = (M 2 — m 2 )/m. 



5 In Eq. ( p)4] ) we used m for the pole mass, but in this section we write m po i , and reserve "m" for a generic 
mass-scheme. 

6 In practice, determining the pole mass from the analysis of experimental data leads to values that depend 
strongly on the order of perturbation theory that has been employed for the theoretical predictions. This 
makes the treatment of theoretical errors difficult. 
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As discussed in Sec. [TIB, it is necessary that the residual mass term is consistent with 



the bHQET power counting, i.e. 

5m ~ §t ~ % ~ T . (98) 

Eq. (|98|) restricts us to a suitable class of short-distance mass schemes for jets. In any short- 
distance mass scheme which violates Eq. (^) the EFT expansion breaks down, and thus 
the notion of a top-quark Breit Wigner distribution becomes invalid. The most prominent 
example for an excluded short- distance mass scheme is the MS mass scheme, m, for which 
?7ipoie — m = 5m. Here 5m ~ 8 GeV 3> T, or parametrically 5m ~ a s m ^> T. Using Eq. ( |96| ) 
and converting to the MS scheme with the 0(a s ) residual mass term we have 

r, r— x if r UsT)5m } 

B +^^ 5m ) = ^( [ ( ^ , p -y + [ ( y , r2]2 |- (") 

Here the first term is ~ l/{mT) and is swamped by the second term ~ a s /T 2 , which is 
supposed to be a perturbative correction. This means that it is not the MS mass that is 
ever directly measured from any reconstruction mass-measurement that uses a top Breit- 
Wigner at some level of the analysis. We stress that this statement applies to any top mass 
determination that relies on the reconstruction of the peak position of an invariant mass 
distribution. 

To define a short distance scheme for jet reconstruction measurements, mj, we choose the 
residual mass term 5mj such that, order-by-order, the jet functions B± have their maximum 
at s t — St — 0, where B + (s) is the gauge invariant function defined in Eq. fl8"8"|). So order- 
by-order in perturbation theory the definition is given by the solution to 

dB + (s, /i, 5mj) 



ds 



= 0. (100) 

s=0 

We call this mass definition the top quark jet-mass, mj(/i) = m vo \ e — 5mj. Since the bHQET 
jet functions have a nonvanishing anomalous dimension, the top jet-mass depends on the 
renormalization scale fi, at which the jet functions are computed perturbatively. Thus the 
jet-mass is a running mass, similar to the MS mass, and different choices for /i > T can in 
principle be made. 

For simplicity we will use the notation B±(s,/j,) for the bHQET jet-functions in the 
jet-mass scheme. At next-to-leading order in a s , 

B ± (s,n) = B ± (s,v) + — { fi r ^ , ( 101 ) 

nmj [s z + 1 A ) A 

where mj = mj(/i) and B + is the pole-mass jet function to 0(a s ). Here we dropped all 
corrections that are power suppressed by T/m. The one-loop relation between the pole and 



jet-mass is |)6] 



mjQi) = m pole - T^P- In (£) + ^ . (102) 



/i\ 3 
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For n = T we have 5mj ~ 0.26 GeV, so the jet-mass is quite close to the one-loop pole mass. 
Equation ( |102| ) also shows that the jet-mass is substantially different from the short- distance 
masses that are employed for tt-threshold analyses [15j , where 5m ~ a 2 s m ~ 2 GeV is of order 
the binding energy of the ti quasi-bound state. Nevertheless, in some of the threshold mass 
schemes |SD| 5m is proportional to a cutoff scale that could in principle be adapted such 
that they are numerically close to the jet-mass we are proposing. A detailed discussion on 
the impact of switching from the pole to the jet-mass scheme at the one-loop level and at 
higher orders will be given in Refs. and [[77], respectively. We remark that many other 



schemes satisfying Eq. fl98|) can in principle be defined, but the existence of one such scheme 
suffices. However, for any suitable short- distance mass scheme the renormalization scale in 
a s contained in 5m has to be equal to the scale \x used for the computation of the bHQET 
jet functions. 

The other function that must be modified in the factorization theorem is 
H m (m, Q /m, fA m , }J,a) ■ However this function only depends logarithmically on m, and 

'»(TH"(t) + °(^)- (103) 

So dropping these perturbatively suppressed power corrections we can simply replace m — > 
mj in H m . We note that any \x dependence from mj(/i) in H m is also power suppressed. 
Thus our final result for the cross-section in terms of the short- distance jet-mass is 



(104) 



x [°°d£ + dr b + (s t - r, fi) b_ (s- t - QL, r , ^ s hcmi (e + , r , ^ , 

where the running jet-mass mj = mj(fi). 
H. Renormalization-Group Evolution 

In order to explain the //-dependence of the factorization theorem in Eq. ( |104j ) we give a 
brief discussion of the renormalization group evolution. A more detailed discussion is given 
in Ref. [BH]. Equation ( |104| ) depends on two renormalization scales, // TO and \i. The matching 
scale fi m ~ m was the endpoint of the evolution of the hard function Hq(Q, fi m ). From the 
matching at m we get the dependence on \i m in H m , and from running below m we get an 



additional dependence on /x as well as Q/m (which is discussed in more detail in Ref. f6~6 
and signifies the presence of a cusp term in the anomalous dimension, see Ref. ||78|| ). The 
/x-dependence in H m cancels against the /i-dependence in the bHQET jet functions and the 
soft function. 

To sum the remaining large logarithms we have in principle two choices. We can either 
run the Wilson coefficient H m , or we can run the individual functions B± and S. The 
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first option essentially corresponds to running the bHQET top pair production current of 
Eq. (|37j), and we will call this method "top-down". The relation 



H m (m, — , fi m , fx) = H m (m, fi m )U Hm (fx m , fx, — ) (105) 

defines the corresponding evolution factor Ug m that is shown in Fig. [5| The second option 
means running the jet functions B± and the soft function Shemi independently with the 
evolution factors Ub ± ((x, fJ> m ) and Us([x, fx m ) respectively, and is also illustrated in Fig. ||. 
This running involves convolutions, such as 

d f 
LM-j-B+{s, = Jds' 7 B+ (s - s') B+(s',fx) , 

B+(s,fj, m ) = Jds U B+ (s - s',n m ,n) B+{s',n) , (106) 

and analogously for B_ and Shemi- Since this method for the running usually involves taking 
the functions B± and Sh em i as an input at the low scale (to avoid the appearance of large 
logs) we will call this option "bottom-up". Because the running of H m is local (i.e. has no 
convolution), this RG evolution only affects the normalization of the cross section and does 
not change the dependence on St and sj in a non-trivial way. This is more difficult to discern 
from the bottom-up running, but when the convolutions for B± and S are combined they 



must become local. These cancellations are discussed in detail in Ref. |66[ where also the 
full leading log evolution is derived. 

Generically, we may wish to run the soft function and jet function to slightly different 
low energy scales. Lets examine the case shown in Fig. |5| where we run the soft function to 
fx a, but run the bHQET jet functions to a slightly lower scale fxr- (The opposite case could 
of course also be realized.) In this case the running is local up to the scale fx a, and below 
this scale we have convolution running for B±. Using Eq. ( |106[ ) the factorization formula 
for split low energy renormalization scales is 

d 2 a TT _ ,__ / Q 



er H Q (Q,fx m )H m (mj, , fx m , i^a) (107) 



dMfdM? ■ v mj 

/oo 
ds' t ds't U B+ (s t -s' t ,fXA,tir) U B _(st-s' f , fi A , A*r) 
-oo 

x / d£Ue-S hc ^(e + J-,fM A )BJs' t -^,T,fj, r ) £_(s;--5L I> r ), 
J-oo \ mj J \ mj J 

where parametrically fx a ~ A*r and here we take mj = mj(/ir). In this paper we will use 
common low energy scales for our numerical analysis, as shown in Eq. ( p.04j ), and leave the 
discussion of the more general case in Eq. (|107|) to Ref. |66| . 
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I. Thrust and Other Event Shape Variables 



Starting from the two-dimensional distribution, d 2 o jdM 2 dM 2 in Eq. ( |104j ) it is straight- 
forward to derive results for other event shape variables. For example, for the thrust T 
defined in Eq. (|59|) , we have 1 — T = (M 2 + M?)/Q 2 which follows using Eq. (gg) with 
Eqs. (|^) and fl60|). Inserting the identity 

M 2 + Ml 



dT6[l-T 1 ^ 2 1 ) (108) 



into Eq. ( |104j) and integrating over M 2 and M 2 we find 

= Voil- 1 ) [ ds t dst B+ (—, T, (i) B_ (—, T, fi) S tbxust (l-T - ( 2m J + s t + s t) ^ ^ 
al j_ OQ \mj / vmj / \ 



(109) 



where <Jq(p) = (ToHq(Q, p, m )H m (mj,Q/mj, p, m , //). Here the thrust soft-function is simply 
a projection of the hemisphere soft function, 

5 thrust (r,/i) = tZT^T - (g+ g £ ^ hcmi^.r,//) (110) 

= ]v-EK r -^^) (0|F ^ n(0)|Xs)(Xs|r " tF - (0)|0) - 

Another well known distribution, which is also frequently analyzed for massless jets, is 
the heavy jet mass. It can be defined by the dimensionless variable 

p = ^Max{M 2 ,M|}. (Ill) 

Using the same steps as above for p, the factorization theorem for top initiated jets is 

— =<7q{ii)1 ds t dst B + (—,T,p) B-l—,T,fj,) Shjm(p - -^,s t ,Si) , (112) 

where the relevant soft-function is 

Sum(p,s h st) =J^d£ + de-5^p-^M^{Q£ + +s t ,Qr+st}^S hemi (e + ,r,fi). (113) 

Factorization theorems for other event shapes that are related to d 2 a/dM 2 dM 2 can be 
derived in an analogous manner. As should be obvious from the definitions of thrust and the 
heavy jet mass distribution in Eqs. (|109|) and ( |112j ), these event shape distributions are also 
characterized by a peak at shape parameter values that are sensitive to the short- distance 
top-quark mass. It is therefore possible to use these event shapes to measure the top-mass 
with a precision comparable to the invariant mass distribution discussed in the previous 
subsection. A brief numerical analysis of the thrust distribution is given in Sec. [IV A. 
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IV. ANALYSIS OF THE INVARIANT MASS DISTRIBUTION 



A. A Simple Leading Order Analysis 

The main result of this paper is the formula in Eq. ( |104| ) for the double invariant mass 
distribution with a short distance top-quark mass suitable for measurements using jets. 
In this section we discuss the implications of Eq. ( |104| ) for top-mass measurements. For 
convenience we rewrite the cross-section in terms of dimension one invariant mass variables 

dPa 4M t Mta? w , , 

F(M t ,Mf, n) , (114) 



dM t dMt (mjT) 2 

where crjf = (ToHq(Q, fi m )H m (mj,Q/mj, fi m , /i) is the cross-section normalization factor 
with radiative corrections, Q is the cm. energy, and we have defined a dimensionless function 

F(M t , M f , n) = {mjYf [°d£ + d£~B + (s t - T, ^)b_ (s- t - 9*-, T, fi)s hemi (i + , £~,fi). 

J-oo \ mj J \ mj J 

(115) 

In terms of M t and Mi the variables s t) i are 

s t = 2M t -2mj, St = 2Mt-2mj, (116) 

up to small Y/m power corrections. In Eqs. (|114[ - |116| ) the jet hemisphere invariant masses 
are M t and Mi and the short-distance top-quark mass that we wish to measure is mj. In 
d 2 a/dM t dMi the function F dominates the spectrum, while 4M f Mf<r^ /(mjT) 2 acts as a 
normalization constant (since M t Mi is essentially constant in the peak region of interest). A 
measurement of the normalization is not optimal for determining mj; it only has logarithmic 
dependence on the short-distance mass, and has larger theoretical uncertainties. On the 
other hand, the spectrum is very sensitive to mj, so henceforth we focus on F{M t ,Mi,fi). 

From Eq. ( |115| ) F is given by the convolution of the computable B± functions, with a 
non-perturbative hemisphere soft-function, S^cmi, ^ na ^ describes soft final-state radiation. 



The majority of the important features of Eq. ( 115 ) can be explained without discussing 
perturbative corrections, so we focus here on the leading order result. From Eq. (|96|) , B± 
are simply Breit-Wigner's at leading order, 

B+(s t ) = ( 1 r s p /T i 2 , : , = / 1 r x p /r ! 2 . ; • ( 117 ) 

ir(mjT) {s t /Ty + 1 TT(mjT) {si/Vy + 1 

For our numerical analysis we use the two-loop standard model prediction for the top-width 
T = 1.43 GeV [7S|| and we take the short distance jet-mass to be fixed at rrij = 172 GeV. As 



demonstrated in Sees. |T| and PH Shemi is the same function that controls the soft radiation 



for massless dijets, which was studied in Refs. p2L 133, p6[. Hence, it is convenient for our 



analysis to adopt the model used to fit the massless dijet data [56 



, M1 ns^^M f£ + £-y-i _ r -(£+y - (£-f - 2b£+t 



cr, n = w(n-^ {-jr ) exp — ^ ) . (118) 
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FIG. 7: Plot of F(Mt, Mf), which is the double differential hemisphere invariant mass cross-section 
d 2 a/dMtdMi in units of A(Tq /T 2 . The observed peak position (intersection of the magenta lines) is 
not given by the true top-quark mass, m = mj = 172 GeV (red lines). This peak shift depends on 
the energy Q, the width T, and the soft-radiation function. The result is shown for Q/mj = 4.33 
and the parameters in Eq. (119). 



Here the normalization constant Af(a,b) is denned so that J d£ + d£~S(£ + ,£~) = 1, the pa- 
rameter A ~ Aqcd sets the scale for £ ± and hence the soft radiation, and the parameter 
a controls how fast the soft-function vanishes at the origin. The dimensionless parameter 
b > — 1 controls the correlation of energy flow into the two hemispheres. Any b ^ implies 
cross-talk between the two hemispheres. A fit to the heavy jet mass distribution using e + e _ 
dijet data from LEP and SLD with Q = uiz gives |5! 



a = 2, b = -0.4, A = 0.55 GeV, (119) 

These values were shown to yield accurate predictions for the heavy jet-mass and C- 
parameter event shapes for a wide range of energies, Q = 35-189 GeV [j5"6|] , as well as 
available thrust distributions with Q = 14-161 GeV [031. We adopt Eq. ( |119| ) as the central 



values for our analysis, but will also discuss how our predictions vary with changes to these 
model parameters. 

In Fig. we plot F(M t ,M { ) using Eqs. (|TT7^TO ) and taking Q = 745 GeV. The key 



feature to note is that the observed peak position is not given by the short-distance top- 
quark mass mj, but is instead shifted upward by ~ 1.5 GeV. The positive sign of this shift 
is a prediction of Eq. ( |115|) irrespective of the choice of parameters. The precise value for 
this shift depends on Q/mj, T, as well as the parameters of the soft function. A less obvious 
feature of Fig. ^ is that the width of the observed peak has also increased beyond the width 
T of Eq. ( |117j ). Physically, the reason for this behavior is that soft radiation contributes to 
the invariant masses, while the Breit-Wigner is only a leading order approximation for the 
spectrum of the top-quark and accompanying collinear gluons. Thus the arguments of B± 
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FIG. 8: Plot of the thrust distribution, da/dT in units of <Jq , for top-initiated events in the peak 



region. We use Q/mj = 5, mj = 172 GeV and the soft function parameters in Eq. (119). 



in Eq. ( |115| ) subtract the dominant soft momentum component from s t f. If we approximate 
<Shemi(^ + , t~) as a very narrow Gaussian centered at = ^q, then the observed peak simply 
occurs at M t) j ~ vrtj + Q£^ /(2mj). Although this model is too naive, we demonstrate in 
the next section that the linear dependence of the peak shift on Q/mj is in fact generic and 
independent of the soft-function parameters. The peak width also increases linearly with 
Q/mj. 

The presence of the shift is due to the inclusion of soft radiation in the definition of 
the invariant masses M t and Mi. Although we adopted a hemisphere mass definition, the 
same type of shift will be present for any jet algorithm that groups all the soft radiation 
into the jets identified for the top and anti-decay products, as we discuss in Sec. [V]. The 
numerical analysis performed in this section applies equally well to these situations, though 
the appropriate definition and model for the soft functions S for such analyses will in general 
be different than that in Eq. ( |118| ) with Eq. ( |119| ). We are not aware of studies where models 
for such soft functions were discussed. 

It is important to emphasize that the shift of the observed peak position away from 
mj is not an artifact of the mass-scheme. At the order used to make Fig. [7| we could set 
mj = m pole since as explained in Sec. [Ill G| they differ by 0(a s T). 7 In a generic short 



distance top-quark jet-mass scheme there is a small shift ~ a s T in the peak position due 
to perturbative corrections in the matrix element defining B± (as discussed in detail in 
Ref. [f}(|). In Sec. |III G| we defined mj using a jet-mass scheme which keeps the peak of B± 



fixed order-by-order in perturbation theory. In this scheme the shift in the peak location 
relative to the short- distance mass is entirely due to the non-perturbative soft radiation. 



7 In general use of TO pole is not a good idea, since in fits it would induce an unphysical change in the required 
parameters a, b, A order-by-order in perturbation theory 
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Although mj is not determined by the peak-position, the shape of the cross-section is 
very sensitive to mj, and hence for precision 5m t < 1 GeV the top-quark mass should be 
determined by a fit to F in Eq. (|114j ). In Sec. [Ill 1| factorization theorems for related event 
shape variables were derived, including thrust da/dT, and the heavy-jet mass da /dp. These 
event shapes also exhibit a peak. They are sensitive to the top-quark mass parameter mj, 
and can be used for top-mass measurements. As an example, in Fig. |8] we plot da/dT 
using Q/mj = 5, mj = 172 GeV, and the parameters in Eq. ( Jl 19 ) . The expected peak 
in the thrust distribution is at 1 — T ~ 2m 2 /Q 2 = 0.08, and is shifted to the right by 
A(l — T) = 1.3 x 10~ 3 by the soft-radiation. A gain the direction of the shift is a prediction, 
but the precise amount of the shift depends on the soft-model parameters in Eq. (|119|) 
as well as Q/mj. An analysis of any other event shape distributions that are related to 
d 2 a/dM 2 dM 2 can be made in a similar fashion. 

In Sec. [IV B| we explore the functional dependence of the peak shift for d 2 a/dM t dMi in 
greater detail. In Sec. |IVC] we discuss the implications of our results for fits to determine 
the short-distance mass. 



B. Analysis of the Peak Shift and Broadening 



In this section we analyze the parameter dependence of the peak shift and broadening of 
the width, and demonstrate that they have a linear dependence on Q. The main analysis is 
carried out assuming that the soft-function model parameters have been determined from 
massless jet observables with small uncertainties and adopting the parameters in Eq. (|119|) . 
It is, however, also instructive to study the dependence of the invariant mass distribution 
on variations of the model parameters, anticipating that the soft function is different when 
the definition of the invariant masses is modified (see the discussion in Sec. |V|). We carry 
out such an analysis near the end of this section. 

In Fig. [|a we plot the peak location, Mf eak , for nine values of Q. M f peak is obtained from 
the two-dimensional distribution, and corresponds to the intersection of the magenta lines 
in Fig. [?|. Since d 2 a/dM t dMi is symmetric the value of Mf eak is the same. Note that for 
Q ~ 2m j where the tops are near ti production threshold, our effective theory expansions 
do not apply. The straight blue line in Fig. |9]a is a linear fit to the points with Q/mj > 4, 
and clearly shows that the peak location grows linearly with Q. In Fig. ||b we plot the "Peak 
Width", defined as the full- width at half- max of d 2 a/dM t dMi in the top- variable M t , while 
fixing the antitop Mj = M? eak . The red solid line is a linear fit for Q/mj > 4. This figure 
demonstrates that we also have linear growth with Q for the width of the measured invariant 
mass distribution. Note that the values for the peak position and peak width shown are 
consistent with our power counting since s t j can be order V as well as greater than V. 

To get a better picture of how the distribution changes with Q we plot the single invariant 
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FIG. 9: Effect of a change in Q on the invariant mass distribution. Results on the left are generated 
from d 2 a / dM t dMi, a) shows the peak position versus Q/mj, and b) gives the full width at half- 
max versus Q/mj. In c) we show da/dM t in units of 2a^ /Y for different values of Q/mj. The 
curves use mj = 172 GeV, Y = 1.4 GeV, and the parameters in Eq. (|ll 



mass distribution da/dM t in Fig. |9]c. In particular we plot 



F\{M t ) 



dMt F(M t ,M { ), 



(120) 



which gives da/dM t in units of 2<Jq /Y. In the numerical analysis we center the integration 
interval [Mi ower , M upper ] on M? eak with a size that is twice the measured peak width. Hence 
the size of the interval depends on Q, but keeps the number of events collected at each Q 
constant for the comparison. For different choices of Q we find that the peak position and 
width of Fi(M t ) behave in an identical manner to Figs. ^|a,b, including having essentially 
the same slopes. In order to keep the area under the curves constant the peak height drops 
as Q is increased. Note that for values Q/mj ~ 8-10 the observed peak location may be as 
much as 2.0-2.5 GeV above the value of the Lagrangian mass mj one wants to measure. In 
our analysis mj is held fixed as shown by the dashed line in Fig. |9|c. 

To gain an analytic understanding of this linear behavior we consider the effect of Q on 
the mean of the cross-section, which is a good approximation to the peak location. Taking 
the first moment with respect to s t /2 = (M t — mj) over an interval of size 2L ^> QA and 
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the zeroth moment in sj/2 = (Mf — mj) gives 

-L 5 /"oo /-oo />L 



= / ds t S -^ [°dstF(M t ,Mt) = rd£ + [ ds t S -±B + (s t -—) [°d£~ S hcmi (£ + , £~) 

mjT 2 J_ L 2J_ 00 J_ L 2 V mjJJ-oo 

- -^-d • (121) 



2m_ 

Thus the mean grows linearly with Q/rrij with a slope determined by the first-moment of 
the soft function, sjH = / £ + 5 hemi (£ + , In the first equality of Eq. the 

.B- function drops out because we integrate over all %. The approximation in Eq. ( |121| ) is 
that terms of ~ 1/L are dropped. We can also directly consider the location of the peak in 
M t , again integrating over Mi for convenience. We use the fact that the tree- level B + (it) is 
symmetric, and solve for M f peak = mj + s peak /2 by setting 

./-oo J -oo ^ m J ' J-oo 



(Sl _^W (0) a( Sl -^f ( o) + ... 

For QA 3> mf we can keep only the first term which yields 



Mr k s£2 



2m, ~ hemi ' 



Aemi^V") • (122) 



(123) 



Thus we find the same shift as for the moment in Eq. (|121|) . Our default model in Eq. (|119|) 
gives S^/2 = 0.31 GeV for the slope in Q/mj. This can be compared with the fit to 
the two-dimensional peak position, Fig. [|a, which gives a slope of 0.26 GeV. The fit to the 
peak position of Fi(M t ) in |^c has a similar slope, 0.25 GeV. Finally, the first moments of 
Fi(M t ) also display linear behavior in Q/mj with a slope of 0.28 GeV. We see that S^^/2 
accounts for the largest portion of these slopes, with the remainder being accounted for by 
other moments. Note that the linear behavior in Q/mj observed in Fig. ^ is much more 
accurate than the statement that S^^/2 determines the proper slope at lowest order. We 
want to point out that a first measurement of the short distance mass could be made using 
a couple of different Q values and a simple extrapolation with Eq. (|123|) . 

Finally we consider the effect on the invariant mass shift from a scan over model param- 
eters. F(M t , Mt) depends on the parameters 

mj, T, /? = -^, a, b. (124) 
mjT 

Here the scale A for the soft-function only shows up along with Q/mj in the effective boost 
parameter (3. To demonstrate that it is (3 that appears in F(M t , Mt), switch integration 
variables to x — £ + /A and y = £~/A, and let s t j = Zt,tT. This yields a soft function 
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FIG. 10: Dependence of the invariant mass distribution on the soft-function model parameters. In 
a) we show 9 models with different a and b parameters, and in b) we show the resulting da/dMt in 
units of 2oq /T. In c) we plot the first moment of the invariant mass distribution, F<^ ,0 \ versus the 
first moment of the soft-function l + , S'[ 1,0 1. In d) we plot the peak position of da/dMt versus 
These plots are made with Q/mj = 5 and mj = 172 GeV. Note that these scans are only relevant 
experimentally if the universality of S'hcmi with massless dijet events is not used to determine the 
parameters. 



A 2 S , homi(Ax, Ay) that is independent of A, and B + (s t — Q£ + /mj) = B + (T(z t — /3x)) which 
is only a function of (z t — f3x) times (Trrij) -1 . Hence F(M t ,Mj) is only a function of j3, 
z t f = {M t f — mj)/r, (mjT), and the model parameters a and b. Hence changing A has 
the same effect as changing Q/mj. Below we will only consider variations of the model 
parameters a and b, while keeping A = 0.55 GeV. 

We generate 9 soft-function models from the intersection of a = {1,2,3} and b = 
{—0.9,0.0,0.9}, and in Fig. |TD|a give the profile of these models by plotting S(£ + ) = 
J d£~ S'hcmi (^ + , Increasing a shifts the distribution to larger average momenta. For 



each model the result for the single invariant mass distribution F±(M t ) is shown in Fig. |i0| b 
with curves of a matching color. We again take the Mi integration interval, [Mi ower , M, 



uppperj ; 
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to be centered on the measured peak, with a size that is twice the measured peak width. 
The peak positions for Fi(M t ) are ordered in the same manner as the peak positions for the 
S(£ + ) models. We note that models with 6 = 0, 0.9 generate smaller peak shifts than those 
with b ~ —0.9. To examine the peak shifts more quantitatively we plot the first moment 
plho] versus the first moment Sj^f^ in Fig. P~D|c. To compute F^ 1,0 ^ we restrict the two inte- 
grals in Eq. ( |121| ) to the same interval choice [Mi ower , M upper ]. From the figure we observe 
that the mean of the invariant mass distributions for different models fall close to a straight 
line. In Fig. [ITJ| d we plot the peak position M t peak for each model versus the first moment 
^hemi! anc ^ observe that the behavior is also quite linear. We conclude that the main effects 
of a, b on the peak shift are controlled by the first moment parameter SwL 



C. Implications for top-quark mass measurements 



In this section we take a step back to consider the more general implications of our method 
for top-quark mass measurements. In a realistic top-quark mass analysis at a hadron collider 
with pp or pp collisions, the set of issues that effect the accuracy for a m r measurement 
and that can potentially be improved by theoretical progress includes: i) the choice of the 
observable to be measured, ii) the top mass definition, iii) hadronization effects, iv) color 
reconnection, v) final state radiation, vi) initial state radiation, vii) underlying events, viii) 
cuts to remove the beam remnant, and ix) parton distribution functions. In our analysis we 
treat e + e~ collisions which allows us to investigate strong interaction effects in categories 
i)-v). We briefly discuss what our result for d 2 a/dM t dMi implies for these uncertainties. 

The main advantage of the factorization approach is that it keeps careful track of how 
changing the observable effects corrections from the other categories. For example, switching 
from invariant mass variables to thrust gives a different function for the non-perturbative 
soft radiation, but the soft-functions in these observables are related by Eq. ( |110| ), and 
one model can be used to fit both of them. In our analysis, the inclusive nature of the 
hemisphere invariant mass observable reduces the uncertainty from hadronization effects. 
In particular it yields jet-functions which sum over hadronic states with invariant mass up 
to ~ mT, and remain perturbatively computable due to the low- momentum cutoff provided 
by the top-width. Final state gluon radiation from the decay products contributes to the 
width, while emission of soft-gluons in the cm. frame organize themselves with radiation 
from the top-quark to give a single universal soft-radiation function. Thus, non-trivial 
color reconnection effects between the decay products are power suppressed. The level of 
control provided by the factorization theorem therefore provides a significant reduction in 
the associated uncertainties. Of course the nature of this control is observable dependent, 
and will undoubtably change in the hadronic collider environment, in particular, hemisphere 
masses are not suitable for mass measurements at the Tevatron or LHC. Nevertheless, we 
expect the control provided by the factorization approach to provide useful applications 
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to this well. Finally, as discussed in detail in Sec. |III G| , the inherent theoretical 

ambiguity in the pole mass, 5m polc ~ Aq CD , can be avoided by switching to a short-distance 
jet-mass raj. This mass definition is suitable for reconstruction measurements in e + e~ 
collisions, and in principle also for pp, and pp collisions. 

For the case studied in detail here, a measurement of d 2 a/ dM t dMi from energetic top-jets 
in e + e~ collisions, there are at least two ways the result in Eq. ( |114| ) can be used to fit for the 
short- distance mass raj. In the first one takes the soft-function model parameters a, b, and A 
from a fit to massless dijet event shapes, and then analyzes d?a/dM t dMi to fit for mj. This 
method makes use of the universality of the soft-hemisphere function Sh ern i(£ + , £~) between 
massive and massless jets. Alternatively, one can vary Q and do a simultaneous fit to mj, 
a, b, and the effective boost parameter (3, to determine the soft-parameters from the same 
data used to determine raj. This may be advantageous if the energy resolution, jet energy 
scale, or other experimental effects have non-trivial interactions with the soft radiation that 
are particular to ti decays. In the next section we consider how the factorization theorem is 
modified by the use of an inclusive kx algorithm rather than using hemisphere masses. 

To conclude this section we note that in our analysis we have not accounted for QED 
effects since in this work we are interested in treating the effects of the strong interactions. 
For a realistic description of experimental data obtained at a future Linear Collider QED 
effects will of course have to be included. Effects from QED can contribute to the categories 
ii), v), vi) as well as to the QED analogues of categories iv) and ix). As such the treatment 
of these QED effects is straightforward and can be included naturally in our factorization 
approach. To be more concrete, the QED analogue to parton distribution functions entails 
to account for initial state radiation, beam strahlung and the beam energy spread through 
a luminosity spectral function which has to be convoluted with the QCD cross section. The 
luminosity spectrum is obtained from analysing Babbha scattering [jSOfl . The effect of soft 
photons showing up in the two hemispheres can be incorporated as additional perturbative 
contributions in the soft function, and the effects of collinear photons can be incorporated 
as additional perturbative contributions to the jet functions in our factorization theorem. 
Finally, the effect of hard photons not alligned with the top and antitop jets is analogous 
to the productions of additional hard jets, which leads to contributions in the hemisphere 
masses away from the resonance region. Compared to the QCD effects treated in this 
work the above QED corrections lead to changes of the invariant mass distribution that are 
suppressed by the small QED coupling. 

V. FACTORIZATION FOR MASSES BASED ON JET ALGORITHMS 

In this work we have defined the top and antitop invariant masses up to now as the 
invariant masses of all particles in the two hemispheres defined through the thrust axis of 
each event, see Fig. |3|. In past experimental studies, on the other hand, a kj- algorithm 
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was employed so that each event results in exactly six jets for the all-hadronic decay mode, 
e + e~ — > ti + X — > 6 jets || |9|. Of these six jets, three jets were combined to the top and 
the other three to the antitop invariant mass. We remind the reader that jet algorithms 
for e + e~ collisions do not need to remove any 'beam remnants', so every final state particle 
of an event is eventually either assigned to the top or the antitop invariant mass. It is 
this crucial aspect of jet algorithms for e + e~ collisions that makes them share a number of 
important properties with the hemisphere invariant masses that we have analyzed so far. 
One of these properties is that having both invariant masses in the peak region close to the 
top quark mass automatically ensures that that the event is dijet-like, such that the EFT 
setup discussed in the previous sections can be applied in the same way as it was for the 
hemisphere masses. 

In this section we show that using a jet algorithm with the property mentioned above for 
the top and antitop invariant mass reconstruction, the double differential top and antitop 
invariant mass distribution in the peak region can be written in the factorized form of 
Eq. ( |104j) , but with a different soft function S(£ + ,£~) which depends on the jet algorithm. 



All other ingredients, the jet functions B± and the matching and evolution factors are 
identical. For the proof we assume that the top and antitop decay jets 8 obtained from the 
jet algorithm can be assigned unambiguously to the top and the antitop, i.e. we neglect the 
combinatorial background. This simplification is possible at leading order in m/Q because 
hard jets from the top decay only have a very small probability of order (m/Q) 2 to show 
up in the hemisphere of the antitop quark, as was already pointed out in Sec. [111E| . The 



analogous statement is of course also true for hard jets from the antitop decay. Moreover 
we assume that the jet algorithm uses simple addition of four-vectors as its recombination 
scheme for merging final state objects. 

The proof can be carried out in the EFT setups described in Sec. [IIIF| and Fig. |5|. 



The crucial point that has to be shown is that, at leading order in the power counting, the 
total n-collinear momentum Px n enters exclusively the top invariant mass, while the total n- 
collinear momentum enters exclusively the antitop invariant mass, just as for the hemisphere 
mass definitions explained in Sec. [Il[D|. Furthermore the prescriptions to determine the 
soft function has to be provided for a given jet algorithm. This corresponds to defining 
appropriate projection operators P a and Pf, in Eq. (f7i|) or equivalently the momenta and 
k h s for each state \X S ). Apart from that, the derivation of the factorization theorem goes 
along the same lines as for the hemisphere case described in detail in the previous sections. 

Concerning the assignment of n- and ra-collinear momenta it is easy to see that the top 
and antitop collinear momenta are attributed correctly to the top and and antitop invariant 
masses since we can neglect combinatorial background for the assignment of top and antitop 



For this discussion we deal with the case that the top quarks decay all-hadronically into jets. However, 
our arguments also work in principle for final states with leptons plus jets. 
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decay jets at leading order. Assuming for example a kr jet algorithm similar to Refs. || U 
where all final state particles are combined to exactly six jets, one can also conclude that at 
leading order in m/Q the n-collinear gluons are properly assigned to the top invariant mass, 
since these gluons are radiated into the n-hemisphere and therefore assigned to one of the 
three hard jets from the top quark decay. The analogous conclusion, of course, also applies to 
the n-collinear gluons. This shows that the top mass reconstruction based on a jet algorithm 
treats n- and n-collinear momenta essentially in the same way as the hemisphere method. It 
also means that the double differential top and antitop invariant mass distribution based on 
a jet algorithm can be derived in complete analogy to the hemisphere case and has the form 
shown in Eq. ( |104|) . The soft function depends on the jet algorithm that is employed, and 
in particular on the distance measure implemented in the algorithm. Whether a soft gluon 
of a given energy ends up contributing to the top or the antitop invariant mass depends 
on its relative angles with respect to the hard jets coming from the top and antitop decay. 
So upon averaging over the hard jet-configurations, a soft gluon with a given energy and a 
given angle with respect to the thrust axis contributes either to the top or to the antitop 
invariant mass governed by a probability function that is determined by the jet algorithm. 
This means that a soft gluon in let's say the n-hemisphere has in general a nonvanishing 
probability to be eventually assigned to the antitop invariant mass. The equivalence of the 
top-down and the bottom-up approaches to the RG evolution in the EFT's used to derive 
the factorization theorem further ensures that the RG running of the soft function for a 
given jet algorithm agrees with the running of the hemisphere soft function Shemh although 
their scale-independent terms differ. (We assume that the jet algorithm is symmetric in 
its treatment of top and antitop final states.) The explicit one-loop expressions for the 
RG running and the scale-independent contributions of the soft function for a general jet 



algorithm and for the hemisphere masses will be given in Ref. |66| 



VI. CONCLUSIONS 

The reconstruction of top quark invariant mass distributions is one of the major methods 
for measuring the top-mass m at present and future collider experiments. Using a sequence 
of effective theories to separate effects at different mass scales we presented an analytic fac- 
torization approach for the top invariant mass distribution in the peak region. To be definite, 
we derived the double differential top/antitop invariant mass distribution d 2 a/dM t dMi in 
e + e~ collisions for cm. energies Q ^> m, where M t j are defined as the total invariant masses 
of all particles in the two hemispheres determined with respect to the event thrust axis. 
The factorization formula is given in Eq. ( |104[) and represents the leading order result in an 



expansion in m/Q and T/m, where T is the top quark total width. 

The factorization formula consists of two jet functions for top and antitop quarks, which 
depend strongly on the top quark Lagrangian mass, and can be computed perturbatively 
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order-by-order in a s . It also involves a nonperturbative soft function that describes the 
momentum distribution of soft final state radiation. Using alternative invariant mass pre- 
scriptions, for which the soft particles are assigned differently to M t and Mj, the same 
factorization formula applies, but with a different soft function. The observable invariant 
mass distribution is obtained from a convolution of the perturbative jet functions with the 
nonperturbative soft function. Through this convolution the energy of the maximum and 
the width of the observed distribution are dependent on the cm. energy Q. For a lowest 
order analysis see Figs. |7| and [], and the accompanying discussion. 

A very important outcome of the derivation is that the soft function for the hemisphere 
mass prescription also governs event shape distributions for massless dijet events for which 
plenty of data has been collected at LEP and previous e + e~ experiments. Since the soft 
function can be determined from these data, it is possible to predict the top invariant 
mass distribution based on the hemisphere prescription as a function of the cm. energy 
Q, the strong coupling a s and the Lagrangian top mass in different mass schemes without 
hadronization uncertainties at leading order in the expansion in m/Q, T/m and AqcdA^- 
In principle, this allows to measure a short-distance top quark mass from reconstruction 
with a precision better than Aqcd- 

We have proposed a new short-distance mass scheme called top quark jet mass which can 
be measured with minimized theoretical uncertainties from data obtained at a future Linear 
Collider and which can be reliably related to other known short-distance masses such as the 
threshold or the MS masses. We also expect that, quite generally, the jet-mass scheme will 
provide an appropriate mass scheme for jet related observables involving heavy quarks. 

The factorization approach developed in this work can be applied to determine the top 
quark mass from reconstruction of the top/antitop quark invariant mass distributions at a 
future e + e~ Linear Collider. The at present most precise method to measure the top quark 
mass at a future Linear Colliner is the threshold scan method. It relies on the determination 
of the hadronic i?-ratio for cm. energies around twice the top mass and will provide a short- 
distance top quark mass measurement with theoretical and experimental uncertainties at the 
level of 100 MeV. Compared to the measurement of the i?-ratio for the threshold scan, the 
reconstruction of the top/antitop invariant mass distribution is without doubt substantially 
more complicated. But it has the advantage that it can be carried out at any cm. energy 
above threshold and that substantially more luminosity can be spent of it. Given that 
our factorization approach allows to control the perturvative and nonperturbative effects 
contributing to the invariant mass distributions we believe that it could eventually become 
a method that competes with the threshold scan. 

The factorization ideas proposed in this work can be applied to mass distributions of 
other final state particles produced in e + e~ collisions in a straightforward manner. Notable 
examples include single top quark production, the production of W bosons or of new heavy 
colored unstable particles such as squarks or gluinos in a certain supersymmetric new physics 
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scenarios. They will also be relevant for predicting invariant mass distributions at hadron 
colliders. However, at hadron colliders there are additional complications that still need to be 
resolved as discussed in Sec. |1V (J| . These include the initial state radiation and incorporating 
parton distribution functions, which lead to a distribution for Q and require modifications 
of the concept of event shapes, the large pr cuts needed to get clear signals away from the 
beam remnant, and the effects of underlying events, which need to be taken into account. 
Finally, the algorithm for defining and measuring the invariant mass of jets that contain the 
top-decay products is different in the LHC environment. We plan to address these issues in 
future work. 
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